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ABSTRACT 


The  Crew  Equipment  Retrieval  System  (CERS)  is  proposed  for  space  station  to 
provide  the  capability  to  rescue  an  EVA  crewman  or  to  retrieve  equipment  inadvertently 
detached  from  the  station.  This  research  is  directed  to  model,  simulate,  and  analyze  at¬ 
titude  control  for  the  Crew  Equipment  Retriever  (CER)  with  and  without  a  target  during 
autonomous  attitude  hold.  Time-optimal  and  weighted-time-fuel  optimal  control  laws 
are  derived  using  Pontrsagin's  .Minimum  Principle.  The  CER  baseline  configuration  is 
analyzed  to  accomplish  some  of  the  attitude  control  trade-off  analyses  planned  for  the 
CER  preliminary  design  phase.  Optimal  thruster  size  and  placement  are  evaluated  for 
three-axis  stabilization.  Control  stability  when  the  moment  of  inertia  tensor  changes 
during  target  capture  is  evaluated  for  several  worst-case  scenarios.  Attitude  control 
performance  results  arc  computed  through  computer  simulation.  Simulation  of  the  CER 
baseline  configuration  shows  it  does  not  provide  effective  control  during  capture  of  a 
worst-case  S50  pound  target. 

A  new  CER  configuration  scheme  is  proposed,  evaluated  and  compared  to  the 
baseline  configuration.  Fuel  optimal  and  end-of-mission  performance  for  the  new  CER 
configuration  is  evaluated.  Simulation  of  the  CER  proposed  configuration  shows  it 
provides  eflectivc  control  during  target  capture  for  modified  locations  in  the  capture  net. 
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THESIS  DISCLAIMER 


The  reader  is  cautioned  that  computer  programs  developed  in  this  research  may  not 
have  been  exercised  for  all  cases  of  interest.  While  ever>'  effort  has  been  made,  within 
the  time  available,  to  ensure  that  the  programs  are  free  of  computational  and  logic  er¬ 
rors,  they  cannot  be  considered  validated.  Any  application  of  these  programs  without 
additional  verification  is  at  the  risk  of  the  user. 
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I.  INTRODUCTION 


A.  THESIS  OBJECTIVES 

The  overall  objective  of  this  thesis  is  to  model,  simulate,  and  analyze  attitude  control 
tor  the  Crew  Equipment  Retriever  or  "CER"  with  and  without  a  target  during  auton¬ 
omous  attitude  hold  for  all  mission  phases  [Ref.  1:  p.  1]. 

Specifically,  the  CER  baseline  configuration  is  analyzed  and  the  following  trade  an¬ 
alyses  are  accomplished; 

•  Optimal  thruster  size  and  placement. 

•  Control  stability  when  the  moment  of  inertia  tensor  changes  for  the  total  system 

during  capture  of  a  target  (maximum  S50  pounds). 

•  Optimal  fuel  performance. 

•  End  of  mission  control  performance. 

B.  CERS  CO.NCEPT  DESCRIPTION 

1.  CERS  Origin  and  Purpose 

In  May  19S7  the  National  .Aeronautics  and  Space  .Adntinistration  (N.AS.A). 
Johnson  Space  Center  Space  Station  Projects  Office  sent  out  a  "Request  for  Proposal" 
(REP)  as  part  of  Space  Station  Work  Package  2. [Ref  2  :  p.  L-2-14a]  The  REP  (including 
an  added  .Amendment  7)  defined  requirements  to  provide  the  capability  to  rescue  an  in¬ 
capacitated  E\'.A  ( external-vehicular  activity)  crewman  and  to  retrieve  equipment 
inadvertantly  detached  from  Space  Station. 

.McDonnell  Douglas  Astronautics  Company  (.MD.AC)  responded  to  the  REP 
with  the  Crew  Equipment  Retrieval  System  or  "CERS".  a  practical,  low  cost  retriever 
concept  in  September  19S7,[Ref  1:  p.  i]  Funding  for  CERS  is  predicted  for  Spring  19S9 
and  the  conclusions  and  recommendations  of  this  research  will  hopefully  assist  Johnson 
Space  Center  and  .MD.AC  in  completing  some  of  the  attitude  control  trade-olT  analyses 
planned  for  the  CERS  preliminary  design  phase. 

2.  CERS  Overall  Mission 

The  CERS  includes  the  following  components: 

1.  A  crew  and  equipment  retriever  (CER)  vehicle. 

2.  .A  retriever  support  station  (RSS). 

3.  ,A  retrieser  berthing  station  (RBS). 


4.  Supporting  Space  Station  capabilities. 

The  major  components  of  the  system  are  shown  in  Figure  1  on  page  3.  [Ref.  1:  p.  l.'^] 
The  Cl-R  is  remotely  commanded  from  inside  Space  Station  by  the  retriever  operator 
or  the  Station  Trafiic  Management  System  depending  on  the  mission  phase. 

The  mission  sequence  of  events  is  shown  in  Figure  2  on  page  4  [Ref.  1  :  p.  S]. 
The  major  mission  events  are  further  subdivided  into  nine  phases; 

1.  Pre-mjssion  Phase  (Storage  and  Standby  Mode) 

2.  Mission  Phase- 1  (Target  Departure  and  Acquisition) 

3.  Mission  Phase-2  (CER  Deployment  and  Rendezvous  Initiation) 

4.  Mission  Phase-3  (Target  Rendezvous  Completion) 

5.  Mission  Phase--  (Target  Capture  and  Station  Rendezvous  Initiation) 

6.  Mission  Phase-5  (Station  Rendezvous  Completion) 

".  Mission  Phase-6  (Terminal  CER  Operations) 

S.  Mission  Phasc-7  (.Airlock  Ingress) 

9.  Post-Mission  Phase  (CER  Refurbishment) 

The  major  mission  phases  are  shown  with  a  time-line  in  Figure  3  on  page  5  and  a  pic¬ 
torial  representation  of  the  overall  mission  profile  is  shown  in  Figure  4  on  page  6.  IRcf. 
1:  pp.  5-6,1"] 

3.  CERS  Capabilities 

CERS  was  proposed  as  a  quick-response  design  and  a  summary  of  its  more  im¬ 
portant  capabilities  are  listed  as  follows;[Ref  1;  p.  9] 

•  Retrieve  an  S50  pound  target  (includes  a  lO-’o  safety  margin),  defined  as  a  lost 
crewmember  or  equipment. 

•  Retrieve  and  deposit  target  into  airlock  within  120  minutes  of  deployment. 

•  Retriever  can  be  activated  and  deployed  without  the  assistance  of  an  EVA 
crewman. 

•  Retriever  senses  own  attitude,  range  and  range  rate  to  target,  and  relays  informa¬ 
tion  to  station. 

•  Retriever  can  be  remotely  controlled  from  any  command  and  control  station  in  the 
Space  Station. 

•  Accomodates  a  worst-case  target  separation  of  3.5  ft  sec. 

•  Retriever  senses  and  controls  its  own  attitude  with  and  without  a  target. 

•  Has  attitude  hold  and  three-a.\is  translation  capability. 


Figure  1.  CERS  INlajor  Comijoiieiits;  From  (Ref.  1  :  p.  15] 


4.  CERS  Major  Systems 

CERS  major  systems  are  described  in  detail  in  Ref.  1  (pp.  20-68).  The  Avionics 
system  is  of  particular  interest  as  it  contains  hardware  and  software  required  to  perform 
attitude  control,  attitude  determination  and  command,  control  and  monitoring  functions 
(among  other  important  functions).  Attitude  reference  determination  for  CERS  is  well 
established  and  detailed  in  Ref  1  (p.  25)  and  will  not  be  discussed  further. 

Attitude  control  of  the  Crew,  Equipment  Retriever  (CER)  is  based  on  a  simple 
switchline  control  system  that  maintains  the  attitude  within  a  specified  deadband  about 
the  conunanded  attitude.  Once  the  target  is  acquired,  the  control  system  uses  a  closed 
loop  television  tracking  processor  along  with  tele-operator  conunands  to  orient  the  CER 
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Figure  2. 


Mission  Sequence  of  Events;  From  [Ref.  1  :  p.  8] 


as  needed  (or  the  spccilk  mission  phase.  The  control  system  approach  is  shown  in  Fig¬ 
ure  5  on  pace  7. 

Velocity  control  reciuircd  Ibr  the  CFiR  to  complete  rendezvous  with  the  tartet, 
to  zero  relative  velocity  when  approaching  close  to  the  target,  and  to  rendezvous  with 
the  Statioit  is  conunanded  by  the  Station  TralTic  Management  System  [Ref.  1:  p.  2S]. 
Discussion  of  the  CUR  velocity  adjustments  required  to  initiate  and  maintain  iranslaiion 
to  acliicve  rendezvous  are  noi  included  in  this  research. 

5.  CER  Baseline  Configuration 

Ihc  initial  CFR  baseline  configuration  proposed  to  provide  the  desired  per- 
Ibrmance  characteristics  during  all  mission  phases  is  represented  by  Figure  6  on  page 
S.  1  his  simple  box  shape  representation  of  the  CER  will  be  used  to  simplify  the  calcu¬ 
lations  describing  its  attitude  dynamics. 

The  characteristics  of  the  baseline  configuration  that  will  be  used  in  the  analysis 
of  the  CER  attitude  control  performance  are  listed  as  follows:  [Ref  1:  p.  24j 

•  i  pounds  total  weight  (combined  weight  of  capture  arms  and  nets  assumed  neg¬ 
ligible  for  computation  purposes). 

•  I  hrcc-axis  (six  degrees  of  fi'cedom)  stabilized  remote  tele-operated  free  flyer. 

•  Uses  24  cold  Nitrogen  gas  (Vj)  jet  thrusters  (two  separate  redundant  systems)  rated 
atl.Uibf 
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Figure  ?.  Major  Mission  Pliases  (nitli  time-Jine);  From  [Ref.  1  :  p,  17] 

•  Attitude  control  accomplislied  by  firing  thrusters  in  pairs  to  produce  couples  about 
any  of  the  three  orthogonal  axes. 

•  Suniniation  of  control  torques  {T  =  r  x  F): 


(  +  -)ROLL  =  3  (ft-lbO 

(  +  -)PITCn  ^7,  =  3  (ft-lbO 

(4--)YAW  ^F,  =  4  (ft-lbO 
•  Moment  of  Inertia  Tensor  [/rorJ  Target): 


(1.3) 


Figure  5.  CER  Attitude  Coniiiiand  and  Control;  From  (Ref.  1  :  p.  27] 


/ i(jj  - 

'I  he  moment  of  inertia  tensor  data  for  the  CER  with  no  target  was  calculated  by  the 
author  assuming  an  850  pound  total  system  weight  synmictrically  distributed  about  the 
center  of  gravity  (weight  ol'the  cajiturc  mechanism  assumed  negligible). 

C.  TMESIS  ORGANIZATION 

In  Chapter  II  the  CER  attitude  dynamics  are  modeled  and  described  by  deriving 
Euler's  moment  equations  (or  the  retriever  with  target  and  representing  them  using  state 
variables.  1  he  moment  of'  inertia  tensors  for  the  CER  with  and  without  a  target  are 
calculated  for  din’erent  worst-case  target  capture  scenarios.  Control  torques  sunuiied 
about  the  shifting  center  of  gravity  are  calculated  for  the  CER  with  and  without  a  target. 

Chapter  111  describes  the  attitude  stabilization  of  the  CER  by  deriving  the  time  op¬ 
timal  and  fuel  optimal  control  laws  that  will  be  used  in  the  model  simulation.  State 
space  plots  of  the  corresponding  switching  curves  are  illustrated,  and  a  deadband  control 
scheme  about  the  origin  is  also  implemented. 

In  Chapter  IV  the  CER  model  is  simulated  for  the  different  worst-case  target  cap¬ 
ture  scenarios.  I'he  CER  baseline  configuration  simulation  results  are  tabulated  and 


4v 

4. 

■  39.6 

0 

0  ■ 

Ar 

= 

0 

55 

0 

4. 

A. 

0 

0 

55 

(in  slug-//^) 


(1.4) 
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analyzed  to  yield  attitude  control  performance.  A  new  CER  configuration  is  proposed 
and  the  new  model  is  simulated.  The  results  are  tabulated  and  compared  to  the  CER 
ba.seline  conliguration.  Fuel  optimal  analysis  of  the  proposed  model  concludes  Chapter 
IV. 

Conclusions  based  on  the  simulation  results  arc  presented  in  Chapter  V,  as  well  as 
recommendations  to  improve  the  design  and  operational  performance  of  the  CER. 
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II.  ATTITUDE  DYNAMICS  FOR  THE  CER 


A.  ATTITUDE  PARAMETERIZATION 

1.  Parameterization  Discussion 

Several  representations  of  the  rigid  body  rotations  about  the  coordinate  axes  for 
the  CER  were  considered,  including  the  use  of  quaternions  as  formulated  by  Hamilton 
in  Ref  3.  The  direction  cosine  matrix  is  well  known  as  a  nine  parameter  coordinate 
transformation  between  the  body  or  spacecraft  axes  and  the  selected  reference  frame 
axes,  and  is  considered  to  be  the  basic  quantity  specifying  the  orientation  of  a  rigid  body 
in  space.  [Ref  4:  pp.  410-420]  Parameterizing  the  direction  cosine  matrix  in  terms  of 
Euler  Symmetric  parameters  incorporates  the  use  of  quaternions,  requires  only  four  pa¬ 
rameters  and  does  not  involve  trigonometric  functions.  However,  the  geometrical  vi¬ 
sualization  of  the  rigid  body  rotations  using  Euler  Synmretric  parameters  is  not  readily 
apparent. 

The  author  decided  to  use  Euler  angle  parameterization  of  the  CER  rigid  body 
rotations.  Euler  angles  are  in  terms  of  three  rotation  angles  about  coordinate  axes,  ro¬ 
tated  in  a  specific  sequence  to  describe  the  rigid  body  orientation  in  space.  They  are  not 
computed  as  easily  as  Euler  Symmetric  parameters  but  are  easy  to  visualize  and  small 
angle  approximations  can  be  used. 

2.  Euler  Angles  for  the  CER 

The  CER  baseline  configuration  shown  in  Figure  6  on  page  8  illustrates  three 
orthogonal  coordinate  axes  with  origin  at  the  CER  center  of  gravity.  Euler  angle  rota¬ 
tion  Id.  I,/,  p)  for  the  CER  is  shown  in  Figure  7  on  page  10  and  the  dynamic  elements 
shown  are  defined  as  follows;  [Ref  4;  pp.  516-523] 

•  6  is  the  roll  angle  (in  radians)  about  X 

•  fj/  is  the  pitch  angle  (in  radians)  about  Y. 

•  is  the  yaw  angle  (in  radians)  about  Z. 

•  L  =  L,  +  L,  +  L.  is  the  angular  momentum  vector. 

•  fj  =  c),  +  (j).  =  d  +  >1/  +  (f)  is  the  angular  velocity  vector  (in  radians  second). 
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Note  ilie  coupling  between  the  angular  momentum  vector  components  as  shown  by  the 
small  vector  groups  at  tlte  end  of  the  coordinate  axes.  Since  the  rate  of  change  of  the 
angular  momentum  vector  is  equal  to  the  sum  of  tlie  applied  torques,  or 


^  Ai 

then  the  sum  of  the  applied  torques  about  each  axis  yields 

T  = 

*  Y 


Ty  = 


(2.1) 


r.= 


A/., 

LM 

L-A.A 

Ai 

Ai 

“T 

At 

AL, 

L,A0 

1  . 

LM 

Ai 

At 

“T 

At 

AL. 

LM 

•  4- 

L,A0 

Ai  Ai 


Ai 


(2.2) 

(2.3) 

(2.4) 
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The  angular  momentum  vector  is  also  defined  as 


T  =  C/]w  (2.5) 

where  tu  is  the  angular  velocity  vector  and  C/]  is  the  moment  of  inertia  tensor  given  by 


I  = 


^xx  ^xr 

lyx  hy  ly: 

L-.X  4  .  4.- 


It  follov.'c  that  the  rate  of  chanec  of  the  angular  momentum  vector  is 


Z.  =  [r]w  (2.6) 

Now  taking  the  limit 

=  I 

/  — +  (1 

for  equations  (2.2),  (2.3),  and  (2.4)  and  also  substituting  equation  (2.6)  yields  Euler's 
moment  equations 

1  r  =  /.vC''  —  uh-f..,  -r  (2.7) 


T..  =  I..C)  —  Q}rL,  +  (Jj.l.r 

»  A  A  »  A 


(2.S) 


T-  ~  /-cj  —  -f  oJx^v  (2.9) 

Combining  terms  and  expanding  a)  and  [/]  while  substituting  equation  (2.5)  into 
equations  (2.7).  (2.8).  and  (2.9)  yields 

^,x.x^\x  ^xv^^y  /j.y  /,,)  (00,00 j,)/„,  (0O,CO,)/y,  +  (OO^OO,)/ 

+  («,..x,)0..  '  '  ■ 

T,  =  Ix-/'>x  +  Av^-’v  +  ■"  ( to .^Oi < -)( /,,  —  l^x)  -  {o)x<-Jx)^x:  ~  +  (Wr‘^v)^.Tv 

-^((O/OJ/,,  (2.11) 


II 


“  Ay)  “  iiOya)y)}^y  -  {u)y(0,)lj^,  +  {O^iOjlry 

Assuming  small  angular  velocities,  the  terms  containing  angular  velocity  products  in 
equations  (2.10).  (2.11),  and  (2.12)  can  be  dropped,  leaving 

Ax  =  4x<^x  +  Ay^y  +  A-AA’  (2. 1 3) 

T.^.  =  l^(Oj^  +  ly^tby  +  (2. 14) 

T,  =  4-  (2. 15) 

Summarizing  equations  (2.13).  (2.14).  and  (2.15)  yields 

"^7-=  [/](“)  (2.16) 


or 


A=c/]-'y  T 


(2.1-) 


liquation  (2.1")  gives  the  angular  acceleration  vector  w  for  the  CER  in  terms  of  its  in¬ 
verse  moment  of  inertia  tensor  values  and  the  sum  of  all  torques  that  act  on  it. 

B.  STATE  VARIABLE  REPRESENTATION 

The  angular  acceleration  vector  components  given  by  equation  (2.17)  can  be  used 
as  control  accelerations  about  their  respective  coordinate  axes,  'fo  elTectively  demon¬ 
strate  this,  equation  (2.17)  is  modeled  in  state  variable  form  as  follows; 

A',  =  e  A',  =  e  A',  =  0  (  =  (iy,) 

Xj  =  \p  Xi  =  <P  AA  ='/'(  =  Wy) 

.V3  =  0  A3  =  0  A3  =  0  ( =  (i>-) 

The  inverse  moment  of  inertia  tensor  is  given  by 


adjC/] 

dct[/] 


(2.  IS) 
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and  using  only  the  matrix  component  subscripts 


:/r’ 


(j.a  cr  -yrrr'l  -  (.yirr  -  .vr-y)  (.vva  r  -  .rrri  ) 

—  p.-.vrr  —  yrr.v)  (.r.vrr  —  .vrr.v)  —  (.v.yir  —  ATcar) 

(r.vr\-  —  lar.v)  —  (.v.rcn  —  .vvr.r')  (.y.vv-v  —  .ri'v.r) 

- = - 1 - i - ^ - LI - :! -  (2.19) 

[  (.v.yrrrr)  —  (.v.yirci  )  -  (.vvy.rrr)  +  (.yrvzj.r)  +  (.r-y.xrv)  —  (.vzia'zjc)  ]]  ^ 


To  simplify  equation  (2.19).  the  following  variables  arc  defined: 


a  =  f.yzz  —  yzrr) 
b  =  (.Yvzz  -  -Y.-ry) 
c  =  (.yiyz  —  .vz;,y) 
d  =  U'.vzr  —  yzz.v) 

t  =  (.V.VZZ  —  .VZZ.Y) 


/=  (.Y.Yyr  —  .Yzy.Y) 
g  =  (^v.Yry-yv.vr) 
h  =  I.Y.Yzy  —  .Yyr.Y) 
/  =  (.v.vrv  — .Yj.y.v) 

/  =  Denominator 


The  matrix  representation  of  the  state  equations  is  given  by 


_ 

■Vi 

J 

.V- 

0 

0 

0 

1 

0 

(! 

0 

0 

0 

0 

1 

0 

■Vj 

0 

(1 

1) 

0 

(1 

1 

.V. 

0 

0 

(1 

() 

0 

0 

0 

(,) 

0 

(1 

0 

(,) 

.V- 

0 

0 

0 

() 

(t 

0 

■V3 

_ 

0  0  0 

0  0  0 

0  ( )  0 

^  -h  _L 

./  J 

—d  C_  —f 

J  J  J 

-h  i_ 

J  j  J 


•  X 

T, 

T, 


The  last  three  stale  equations  are  meaningful  and  arc  given  as 


( J  -  '  A'i  =  di  7-,.  _  A  7;  +  i_ 

J  "  J  ■  J  ‘ 


(2.2n) 


(^  =  j  ,v,  =  -  dl  r,  +  ^  t;  -  -  7, 

J  J  ■  J  ' 


(2.21) 


(o  =  )  X\ 


—  T. - -  7,  + 

./  J 


(2.::) 


Equations  (2.20),  (1.21).  and  (2.22)  represent  the  CER  angular  acceleration  compo¬ 
nents  for  three-axis  attitude  control.  The  following  sections  will  yield  the  moment  of 
inertia  tensors  and  control  torque  values  needed  in  the  equations. 
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C.  MOMENT  OF  INERTIA  TENSOR  ANALYSIS 
1.  CER  Without  Target  Moments  of  Inertia 

The  CER  moment  of  inertia  tensor  (without  a  target),  was  listed  in  equation 
(1.4)  on  page  7  and  is  repeated  here 


rw:  = 


39.6  0  0 

0  55  0 

0  0  55 


(slug-//^) 


(2.23) 


Notice  the  cross-products  of  inertia  are  zero  due  to  the  symmetrical  distribution  of  mass 
around  the  CER  coordinate  axes.  The  moment  of  inertia  tensor  was  approximated  by 
the  author  usina 


I  2 1  — c 

I  I  -  ) 


=  >  .U/Tr(P^7;.'')l-;^)^') 


(2.24) 


where  Tri.vi  =  Trace  of  .v.  "I"  is  the  identity  matrix,  and 


p  cix  dy  d: 


(2.25) 


where  p  is  the  density  in  slugs  JY.  (Ref  5:  pp.  236-240]  As  an  example.  /„  is  calculated 
using 


where 


P  = 


mass 

vol 


2 

-I-  r  )p  dx  dy  dz 


S.y.)  32.2 
36 


slugs 


resulting  in 

/„  =  39.6  slug-//' 


(2.26) 


(2.2i 


2.  Target  Moments  of  Inertia 

The  CER  will  be  used  to  capture  various  targets  with  difTerent  mass  and  mo¬ 
ment  of  inertia  properties.  The  flexible  capture  net  will  provide  a  compliant  capture 
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wliile  niaiiuaiiiing  sufTicicnt  rigidity  to  secure  the  target  and  maintain  its  center  of  gravity 
location  [Ref.  1:  p.  53]. 

A  point  mass  \vith  no  moment  of  inertia  properties  is  used  to  represent  targets 
such  as  small  tools  or  equipment.  The  primary  target  of  interest,  however,  is  an  EVA 
crewntan  which  includes  the  combined  mass  of  the  man,  his  environment  suit  or  EMU 
{extra-vehicular  mobility  unit),  and  the  MMU  (manned  manuevering  unit).  [Ref.  6]  For 
a  worst-case  target  capture,  an  EV.\  crewman  would  have  the  following  approximate 
mass  and  moment  of  inertia  properties: 

226  100  percentile  male 

25S  EMU 

33S  M M U  (full  propellant) 
lo  Ancillary  equipment 
12  1.5"..  safet)  margin 

S5u  pounds 

'  27.21  0  0  1 

C I  rRt'ooRY  1  =  0  44.78  0  (slug-//*)  (2.28) 

fARUEr  |_  0  0  41.62  _ 

'I  he  moment  of  inertia  tensor  calculated  in  equation  (2.28)  assumes  the  target 
axes  would  be  aligned  with  the  CER  coordinate  axes  during  target  capture.  For  a  more 
realistic  analysis,  the  target  is  rotated  about  its  center  of  gravity  to  give  a  new  moment 
ol'  inertia  tensor  which  is  found  by  using  the  coordinate  transformation  theory  detailed 
in  Ref  4  (pp.  "01 -“65).  To  summarize 

r'  =  CT>  (2.29) 

=  (2.30) 

where 

r  =  cartesian  vector  for  target 

C  rj  =  transformation  matrix  or  direction  cosine  matrix 

/  '  =  cartesian  vector  for  rotated  target 

and  from  equations  (2.24)  and  (2.25) 
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c/] 


(2.31) 


Multiplying  both  sides  of  equation  (2.31)  by  and 

!:7-][;/][:7]^= j(C7'](Tr(;^;^^))im^-  (2.32) 

moving  the  T's  inside  the  trace  and  substituting  from  equations  (2.29)  and  (2.30)  yields 


(Tr)/-,.'/;'^)!  - 


and  since 


TrinW^)  = 


then 


crjcncTj^ 


I':  >'! 


^)din 


or 


Liy  =  LnLnuy  (2.331 

From  equation  (2.33).  the  original  moment  of  inertia  tensor  for  the  target  is  multiplied 
by  the  transformation  matrix  ^direction  cosine  matrix)  and  its  transpose  to  yield  the  ro¬ 
tated  target's  tensor.  The  general  direction  cosine  matrix  for  a  3-1-2  Euler  angle  rotation 
is  given  by 

ary4>.eM  = 

cos  i)/  cos  —  sin  (1  sin  i/^  sin 
—  cos  6  sin  <? 

sin  il/  cos  d  -I-  sin  0  cos  i//  sin  4> 
where 


cos  <'/  sin  4>  +  s'mB  sin  i//  cos  (f> 
cos  0  cos  (j) 

sin  i/^  sin  </)  —  sin  6  cos  i//  cos  (/> 


—  cos  0 


(2.34) 
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3  -  (Yaw)  rotation  by  </>  about  f 

1  —  (Roll)  rotation  by  9  about  x' 

2  —  (Pitch)  rotation  by  ij/  about  r" 

Using  equations  (2.33)  and  (2.34)  lor  a  3-1-2  rotation  angle  sequence  of 

0  =  45° 
e  =  A5° 

0  =  45° 


vields 


[/]' 


43.6 

2.4 

-1.4 

2.4 

38. S 

6.4 

-1.4 

6.4 

31.2 

(2.35) 


Out  of  several  dilTerent  combinations  of  rotation  angles  calculated,  equation  (2.35)  re¬ 
presents  the  worst-case  moment  of  inertia  tensor  for  a  rotated  primary  target. 

3.  CER  With  Target  Moments  of  Inertia 

The  CER  with  target  moment  of  inertia  tensor  is  calculated  using  the  parallel- 
axis  theorem 

C/7or]  =  +  C/r..;tc7£r3 

where  is  the  skc'.v-symniciric  niairix  derived  from  the  vector  /s  between  the  center 
of  gravities  of  the  CER  (mass  .)/,)  and  the  target  (mass  .lA)  and  is  illustrated  in 
Figure  S  on  page  IS.  [Ref  7  j 

During  target  capture  the  total  moment  of  inertia  tensor  changes  accordingly, 
depending  on  the  target's  moments  of  inertia  and  location  in  the  capture  net.  To  provide 
the  worst-case  (largest),  moment  of  inertia  tensors  for  input  into  equations  (2.20).  (2.21). 
and  (2.22).  the  following  cases  were  analyzed; 

•  C.ASE  \A  (-V^3,)  -  850  pound  point  mass  captured  on  net  edge  located  at 

=  9(/)  +  OiJ)  A-  1.5(A)  from  CER  center  of  gravity. 

•  C,‘\SE  IB  (F„„)  -  S50  pound  point  mass  captured  on  net  edge  located  at 
F  =  5.5S3(/)  -I-  1.5(^/)  +  1.5(A)  from  CER  center  of  gravity. 

•  C.\SE  1C  (Z„,3,)  -  S50  pound  point  mass  captured  on  net  edge  located  at 

/■.  =  5.5S3(/|  ()iy)  -f-  5(A)  from  CER  center  of  gra\ity. 
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|— (VECTOR)  r2  =  A(i)  + 

B(i)  +  C(k) 

L  [R2]  = 

0  -C  B 

C  0  -A 

-B  A  0 

Figure  8.  Skew-Symmetric  Matrix 


AAA  AAA 

•  CASH  1 D  (.V  =  Z.  >’  =  1 )  -  850  pound  point  mass  captured  at  X  =  Z.Y  —  1  from  net 
center,  located  at  =  7.097(0  -I-  I(/)  +  3.014(A:)  from  CER  center  of  gravity. 

•  CASE  2  (man  &  nunu)  -  850  pound  primary  jarget  (man  &  mmu  rotated 

(f)  0  =  tj/  =  45°)  captured  at  net  center,  located  at  =  5.583(0  +  0(/)  -I-  I.5(A:)  from 

CER  center  of  gravity. 

The  target  capture  cases  are  illustrated  in  Figure  9  on  page  19. 

An  example  calculation  of  the  total  moment  of  inertia  tensor  for  CASE  2  (man 
lC  mmu)  is  in  order.  The  following  values  are  used  in  equation  (2.36); 

Tj  =  5.583(0  +  0(J)  +  l.5{k)  ^ 


where 


I 
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Figure  9,  Target  Capture  Case  Locations 


A/,  .1/2  (26.3975)(26.3975) 

A/,  +  A/j  "  (26.3975)  +  (26.3975) 

=  equation  (2.23) 

Utarcet^  =  equation  (2.35) 


=  13.19S75  (slugs) 


to  yield 


■  39.6 

0 

0  ■ 

'  43.6 

2.4 

-1.4  ■ 

■  -29.7 

0 

110.5 

0 

55 

0 

-f 

2.4 

38.8 

6.4 

- 

0 

-441.1 

0 

0 

0 

55  _ 

-1.4 

6.4 

31.2 

110.5 

0 

-411.4 

or 


112.9  2.4  -111.9 

2.4  534.9  6.4 

-111.9  6.4  497.6 


All  of  the  moment  of  inertia  tensors  for  the  CER  baseline  configuration  during 
target  capture  were  calculated  as  in  the  previous  e.xample  and  are  shown  in  Figure  10 
on  page  21.  A  quick  comparison  shows  that  the  850  pound  point  mass  (with  no  huernal 
cross-products  of  inertia),  induces  larger  cross-products  of  inertia  when  captured  at  the 
edges  of  the  net  than  does  the  primary  target  (850  pound  man  &  mmu),  when  captured 
in  a  rotated  position  at  the  center  of  the  net.  The  ability  to  reduce  the  cross-products 
of  inertia  for  the  CER  (with  target  )  will  prove  to  be  a  critical  factor  for  attitude  control 
stability. 

D.  CONTROL  TORQUE  ANALYSIS 

The  CER  baseline  configuration  was  listed  in  Section  1.5  on  page  4  as  having  24  cold 
gas  jet  thrusters  rated  at  1.0  Ibf  to  make  up  two  separate  redundant  systems.  Also  listed 
were  the  summation  of  control  torques  based  on  the  thruster  size  and  locations  from  the 
CER  center  of  gravity.  The  thrusters  are  fired  in  pairs  to  produce  couples  about  any 
of  the  three  orthogonal  axes.  An  example  of  a  thruster  pair  firing  to  produce  a  positive 
roll  about  X  is  shown  in  Figure  1 1  on  page  22.  The  corresponding  calculations  of  the 
control  torques  summed  about  the  center  of  gravity  are  as  follows: 


F,  =  0(0  +  0(/)  +  \{k) 


(2.37) 
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Figure  10.  CER  Baseline  Configurafioii  Moment  of  Inertia  Tensors  for  Target 


Capture 


=  -2{/J+  1.5(/c) 


(2.38) 


T2  =  0(/)  +  0(/)- l(/<)  (2.39) 

r,^-2(i)-\.5ij)+l.5(k)  (2.40) 
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Substituting  equations  (2.37)  thru  (2.40)  into 

7i  =  qx/)  (2.41) 

and  summing  the  torques 

V  7  =  r,  +  =  (1.5(/)  +  2(/-)  +  0(A))  +  (1.5(/)  -  2(j)  +  0(A))  (2.42) 

or 

^7-^  =  3(0(ft-lbO  (2.43) 

Equations  (2.17),  (2.43).  and  the  moment  of  inertia  I„  value  from  equation  (2.23)  for  the 
CER  with  no  uirgct  are  used  as  follows  to  find  the  normalized  control  acceleration  if  no 
cross-products  of  inertia  exist: 
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(ft-lbO 

(slug-//^) 


(4). 


39.6 


= .07575 


rad 


sec 


(2.44) 


Equation  (2.44)  reveals  the  relatively  small  magnitude  of  the  control  acceleration  gener¬ 
ated  by  the  1.0  Ibf  baseline  thrusters  acting  on  the  CER  with  no  target. 

During  target  capture  the  moment  of  inertia  tensor  changes  and  the  center  of  gravity 
shifts  accordingly.  Since  the  thrusters  are  fired  in  pairs  to  produce  couples  about  any 
of  the  orthogonal  axes,  ihe  shifi  in  center  of  gravity  has  no  effect  on  the  summation  of 
control  torque  calculations.  To  change  the  control  torques,  either  the  location  of  the 
thrusters  must  change  to  increase  or  decrease  the  distance  between  the  thruster  pairs, 
or  the  thruster  size  must  be  increased  or  decreased  as  desired.  The  ability  to  change  the 
control  torques  (and  subsequently  the  control  accelerations),  will  prove  to  be  beneficial 
to  the  attitude  control  stabilitv  of  the  CER  with  and  without  a  taraet. 


III.  ATTITUDE  STABILIZATION 


The  CER  will  provide  an  autonomous  attitude  hold  capability  at  all  times  and  will 
also  respond  to  commanded  rotations  as  required  for  the  mission.  Since  the  CER  will 
operate  in  close  proximity  to  the  Space  Station  and  must  be  precisely  oriented  during 
target  capture,  a  minimum-time  control  response  is  desired  to  maintain  strict  attitude 
control  during  all  mission  phases.  However,  maintaining  strict  attitude  control  while 
coasting  between  orbits  is  not  critical,  especially  without  a  target,  and  a  minimum-fuel 
control  response  could  be  used.  The  following  sections  will  derive  the  optimal  control 
laws  to  obtain  a  minimum-time  and  a  minimum-fuel  control  response.  The  optimal 
control  laws  will  be  used  in  combination  to  define  switching  curves  for  a  weighted-time- 
fucl  optimal  control  response  that  can  be  used  to  support  a  specific  CER  mission  phase. 

A.  TIME  OPTI.MAL  CO.NTROL  LAW 

The  objective  of  minimum-time  control  is  to  transfer  a  system  from  an  arbitrarx'  in¬ 
itial  state  to  a  specified  final  state  as  quickly  as  possible.  The  optimal  control  system  to 
achieve  a  minimum-time  response  is  normal.y  referred  to  as  a  bang-bang  system,  wliere 
the  optimal  control  switches  between  its  maximum  and  minimum  magnitudes.  (Ref  S: 
p.  25^1 

The  control  acceleration  for  each  axis  will  be  calculated,  as  in  equation  (2.44).  for 
the  CER.  Each  acceleration  (eO  will  be  used  to  drive  its  respective  angular  position  and 
velocity  from  arbitrary'  initial  conditions  to  final  states  of  zero.  An  example  problem 
formulated  by  Kirk  in  Ref  S  :  pp.  249-254,  of  deriving  the  optimal  control  law  and  cor¬ 
responding  switching  curves  can  be  used  for  the  CER  system  to  find  the  optimal  control 
to  accomplish  this  in  minimum-time.  To  summarize  the  derivation  to  follow,  the 
underlying  optimal  control  principles  must  first  be  presented. 

The  optimal  control  (u‘).  is  defined  as  the  control  acceleration  that  minimizes  the 
designated  performance  function  The  minimum  time  performance  function  is: 

n. 

4)=  dt  =  tj-k,  (3.1) 
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where  it  is  assumed  that  the  final  time  is  free  (not  specified),  and  the  final  state  is  fixed 
(specified  at  zero).  These  boundan,'  conditions  are  used  to  help  specify  the  necessary 
conditions  for  optimal  control. 

The  necessary  conditions  consist  of  a  set  of  2n,  first-order  differential  equations  re¬ 
ferred  to  as  the  state  and  costate  equations,  and  a  set  of  algebraic  relations  that  must 
be  satisfied  during  the  control  interval.  The  costate  is  described  using  lagrange 
multipliers  and  the  state  x‘{t),  is  from  the  state  variable  representation  of  the  system. 
The  solution  of  the  state  and  costate  equations  contain  2n  constants  of  integration  that 
are  evaluated  by  using  (2n  +  1)  equations  since  the  final  time  r^is  not  specified  [Ref  8: 
p.  200].  It  will  be  convenient  to  use  the  function  Jf,  called  the  Hamiltonian,  defined  as 

■jf  {x{t).u[t).p{i){i)j)  =  g(x{i),u{i)j)  -T  p\i)\ia{x{i)M(t),i)']  (3.2) 

The  two  variables  in  equation  (3.2)  not  yet  described  are  g,  which  is  the  integrand  of  the 
performance  function,  and  a.  which  represents  the  state  equations.  The  necessary  con¬ 
ditions,  for  all  t  6  ftn.oi],  tire  [Ref  S  :  p.  188] 

(.V  {t).u{t),p  (t)j)  (3.3) 

p  (t)  =  -  (x'(t)M\t).p'(t)j)  (3.4) 

0  =  (.v’(/).«*(0.rV).0  (3.5) 

and 

(.V  (tjltj)  -  p  ((;-)  ]^J^'[x'{tj).u  {tf),p'' (tj),tj)  +  (/(ty),r/-)]<)t/=  0  (3.6) 

Since  the  final  time  (t^)  is  arbitrary  and  the  final  state  is  fixed  at  zero,  then  dXf=  0  and 
equation  (3.6)  yields 

(tj).p\tf),t^)  +  (x‘(ty),ty)  =  0  (3.7) 

In  addition  to  the  necessary  conditions  discussed  above,  an  optimal  control  must 
satisfy  Pontryagin's  minimwn  principle  which  states  that  an  optimal  control  must  mini¬ 
mize  the  Hamiltonian.  [Ref  8:  pp.  227-234]  The  minimum  principle  is  appropriate  to  use 


when  the  admissible  controls  are  constrained  by  certain  maximum  values  (in  this  case  the 
thrusters  are  limited  by  their  maximum  force  settings). 

Applying  Pontryagin's  minimum  principle  modifies  equation  (3.5)  to  yield 

yfi-x  (t).u  {i).p  {i).t)  <  .??■(-'•'  {t),i)  (for  all  admissible  u(t))  (3.S) 

where  u'(t)  is  a  control  that  causes  .>i^(.r'(/),i/(0,pV)d)  to  assume  its  absolute  minimum. 

To  find  the  optimal  control  law  and  corresponding  switching  curves  for  minimum¬ 
time  control  about  one  axis  of  the  CER,  the  system  is  defined  by  the  state  equations 

jf,(/)  =  .r^(t) 

(3  9) 

which  will  be  transferred  from  initial  conditions  to  final  states  at  the  origin  in 
minimum-time  by  the  optimal  control  constrained  by 

I  ?/,!/)  1  <  ;/;(/)  (3.10) 

where  «^(/)  is  the  control  acceleration  previously  calculated  in  equation  (2.44).  and  t/Jt) 
is  defined  as  an  admissible  control.  The  system  is  assumed  to  be  completely  controllable 
and  normal  (no  singular  intervals  exist).  The  applicable  necessar>'  conditions  to  be  used 
are  Irom  equations  (3.3).  (3.4),  (3.5),  (3.S),  and  (3.7).  The  performance  function  to  be 
used  is  from  equation  (3.1 ). 

The  Hamiltonian  for  the  CER  system  is 

J/'tXit).n(tl.pltll  =  1  -f  Pi(t)X2(f)  +  P2{0^a(^)  (jTI  ) 


and  the  minimum  principle  produces 

1  -rp-(t)X2(t)+P2(t)U  (t)<  1  +p^/).V2(/) -H/’2(0Wa(0 


(3.12) 


or 


or  from  equations  (3.13)  and  (3.10)  the  optimal  control  must  be 


n(t) 


—  ujj),  foTp2{i)>0 
for/’2(0<0 


(3.13) 


(3.14) 
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Using  the  Hantiltonian  in  the  necessan.'  condition,  equation  (3.4).  yields  the  costate 
equations 


/'i  (^) 


CX-. 


(3.15) 


/’2(0  =  - 


Ml 

CXi 


=  ~Pdi) 


(3.16) 


The  costate  equation  solutions  are  found  by  integrating  both  sides  of  equations  (3.15) 
and  (3.16)  to  yield 


/>!(0  =  Ci 

(3.17) 

/>;(/)  =  -  Cj/  -y  Co 

(3. IS) 

Equation  (3.18)  indicates  iliat  and  therefore  from  equation  (3.14).  u{r)  can  change 
signs  at  most  of::  time. 

Segments  of  the  optimal  trajectories  can  be  found  by  integrating  the  state  equations 
defined  in  equation  (3.9),  and  after  a  few  manipulations  the  result  is 

(t/  =  +  t/,(M)  (3.19) 

-V-'l  /) 

xAi)  =  -~ - F  (u  =  -u^(i))  (3.201 

Equations  (3.19)  and  (3.20)  each  define  a  family  of  parabolas  and  are  shown  on  a  state 
space  plot  in  ITgure  12  on  page  28  with  the  arrows  indicating  direction  as  time  increases. 
[Ref  S:  pp.  252-253]  By  setting  c<  and  equal  to  zero  in  equations  (3.19)  and  (3.20),  to 
satisfy  the  final  time  condition  on  the  state,  the  state  space  plot  shows  one  continuous 
curve  through  the  origin  as  shown  in  Figure  13  on  page  29. 

Since  one  switching  can  occur  at  most,  the  optimal  control  found  in  equation  (3.14) 
can  be  one  of  the  following; 


Figure  12.  Trajectories  for  U  =  +  Uc  and  U  =  -Uc;  From  [Ref.  8  :  p.  252) 

(Form  1)  +  I  e  [rg.f  ],  or 
.  (Form  2)  -  f  e  or 

w  (/)  =  .  (3.21) 

(Form  3)  +  i  e  and  —  n^(i),  t  e  [r,,r  ),  or 

(Form  4)  -  t/<.(/).  te[/o,/|)and  +«<,(/),  re  [/,,/) 

Direct  minimization  of  tlie  performance  function  yields  the  optimal  controls  for  dilTerent 
initial  conditions 

•  Form  1  ”  Initial  states  lie  on  segment  A  -*  O  with  u'(0  =  • 

•  Form  2  --  Initial  states  lie  on  segment  D  0  with  u'{i)  =  — . 
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Figure  13.  Trajectories  for  U  =  (  + -)Uc  through  the  origin;  From  (Ref,  8  :  p,  253] 


•  Form  3  --  Initial  states  lying  below  both  segments  A  ->■  O  and  B -*  O  will  have 
;/■(/)  =  +1/,  until  Zi  ->  O  is  reached,  then  u\t)  =  — thereafter. 

•  Form  4  --  Initial  states  lying  above  both  segments  /I -♦  O  and  Z? -♦  0  will  have 
!/(/)  =  -It,  until  A  0  reached,  then  it'{t)  =  +u,  thereafter. 

Optimal  trajectories  for  several  initial  state  values  are  shown  in  Figure  14  on  page  30. 
1  he  continuous  curve  represents  the  locus  of  points  where  control  switches  from 
H —  tij!)  to  — h  !t,{/)  and  is  known  as  a  sw  itching  curve.  The  equation  describing  the 
curve  also  defines  the  optimal  control  law  for  minimum-time  control  and  is  given  by 


.Y2(/)|.V;(/) 


(3.22) 


To  simplify  the  use  of  equation  (3.22),  a  switching  function  s(.x(t))  is  defined 


5(.\(t))  =  .r,(0  -P 


■T2(/)!.Y2(/)| 

2u,(0 


(3.23) 


and  implies  that  for 

•  v(.\(t))  >  U  .v(r)  lies  above  A-O-B. 

•  v(.\(f,))  <  0  .r(t)  lies  below  A-O-B. 

•  s(.\(l|)  =  0  .v(/)  lies  on  A-O-B. 
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Figure  14.  Optimal  Trajectories  for  Different  Initial  State  Values;  From  [Ref.  8  ; 
p.  253] 

Tlie  switching  function  is  used  in  the  simulation  program  discussed  in  the  ne.\t  chapter. 


B.  WEIGHTED-TIME-FUEL  OPTIMAL  CONTROL  LAW 

Tlie  CHR  may  have  an  opportunity  to  relax  its  control  response  time  while  coasting 
between  orbits,  consequently  conserving  fuel.  Minimum-fuel  control  response,  however, 
is  not  desired  since  the  final  time  i,  is  not  specified  and  the  system  would  virtually  take 
an  infinite  amount  of  time  to  reach  the  desired  final  state  at  the  origin.  Using  a 
wcighicd-iime-fuel  performance  function  would  allow^  the  CER  to  use  a  minimum-time 
control  response  or  to  conserve  fuel  while  maintaining  a  satisfactory  time  response  as 
required  for  the  mission  phase. 

Tlie  same  method  that  described  the  minimum-time  control  law  will  be  used  to  find 
the  optimal  control  law  and  corresponding  switching  curves  for  a  weighted-time-fuel 
control  response.  The  state  equations  for  the  CER  defined  in  equation  (3.9)  will  be  used 
and  the  control  is  again  constrained  as  in  equation  (3.10).  The  final  time  is  not  spec¬ 
ified  and  the  system  is  to  be  transl'erred  to  a  final  state  x[i^  =  0  while  minimizing  the 
following  performance  function: 


J(U]=  ["[l  +/|!/,(nl]^/r  (3.24) 

Note  that  for  /.  =  0,  equation  (3.24)  becomes  the  minimum-time  performance  function 
as  in  equation  (3.1).  For  /  =  oo.  equation  (3.24)  becomes  a  minimum-fuel  performance 
function  corresponding  to  an  infinite  time  to  reach  the  final  state. 


Beginning  the  derivation,  the  Hamiltonian  is 

=14-/1  u^{t)  I  +  p-Xt)X2{l)  4-  P2{t)uJ,t)  (3.25) 

The  costate  equations  are 

/>;(/)  =  0  (3.26) 

pl(i)  =  -p'(i)  (3.27) 

and  equations  (3.26)  and  (3.27)  have  solutions  of 

p.[i)  =  c.  (3.2S) 

/>:(/)=-  C;/  4-  C;  (3.29) 

Pontryagin's  minimum  principle  yields 

/  I  u  {!)  I  -f  p^dji^  I/)  <  /  I  I  4-  P2(0^a(0  (3.30) 

or 

(/  -!-  P2(n)l  u  (/)|  <  |z/,(f)l(/  -|-P;(/))  (3.31) 

From  equations  (3.30)  and  (3.10)  the  optimal  control  must  be 

-f  ujr)  for  /7;(r)  <  —  / 

0  for  —  /  <  p2(0  <  )■ 

w  (/)  =  -  u^{i)  for  /  <  ^2^0  (3.32) 

undetermined,  but  >  0  for  p2{i)  =  —  / 

flr 

undetermined,  but  <  0  for  f>2^r)  =  +  / 
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Equations  (3.2S)  and  (3.29)  for  the  costate  solutions  show  that  pj(/)  can  change  sign 
at  most  once,  so  the  optimal  control  must  be  one  of  the  following: 


0 

+ 

-  «c(f) 

(0,  +  u^(r)) 

(0,  - 

( +  «c(^).0) 

( -  i<AO-0) 

(  +  u,(r).0.-uAO) 
(  -  m(r),0,  +  uAO) 


(3.33) 


The  control  intervals  that  end  with  u'(/)  =  0  cannot  be  optimal  since  the  final  condition 
cannot  be  achieved  without  the  application  of  a  control.  Therefore 


u'(0 


+ 

- 

(0,  +  u^(r)) 

(0,  -  uAO) 

(  +  uA/).0,  -  w,(/)) 
(  -  uAO-0,  +  w^(rj) 


(3.34) 


and  for  the  controls 


«*(/) 


- 

(0.  -  «,ff)) 

(  +  w,(0.0,  -  uAO) 


(3.35) 


the  trajectories  resulting  must  terminate  at  the  origin  with  an  interval  of  u’(i)  =  —  u^(i) 
control,  therefore,  the  terminal  segments  of  these  trajectories  all  lie  on  curve  5  -+  O  as 
shown  in  Figure  14  on  page  30.  When  «^(/)  =  0,  the  state  equations  become 

-v,(0  =  ->^2W  (3.36) 

X2(r)  =  0  (3.37) 


which  implies 
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.xylr)  =  t'3 

(3.3S) 

.r,(/)  =  C3/  +  C4 

(3.39) 

Time  increases  and  jc,(f)  increases  or  decreases  depending  on  if  is  greater  than  or  less 
than  zero  when  control  switches  to  u^i)  =  0,  and  the  trajectories  are  shown  in 
Figure  15  on  page  34. 

The  trajectories  for  u'(r)  =  (  +  —  )u^(r)  are  the  same  as  shown  before  in  Figure  13  on 
page  29.  The  two  types  of  trajectories  are  combined  in  Figure  16  on  page  35.  Noting 
the  labels  in  Figure  16  on  page  35,  the  optimal  control  laws  are  derived  as  follows: 

=  time  optimal  control  switches  from  +  1  to  0  (somewhere  on  C-K) 
h  =  time  optimal  control  switches  from  0  to  -1  (somewhere  on  K-0) 


On  segment  K-0 


and  integrating 


.Vj(/)  =  x.ir) 


(3.40) 


(3.41) 


yields 

.r*(/2)  =  .v*f/,)  -f.v;(;,)[f2  -  /,] 


(3.42) 


Equations  (3.29)  and  (3.32)  indicate 

/)2(/l)  =  -  c,/,  +C2  =  -/: 

(3.43) 

Pi{h)  =  -<^\h  +  C2  =  +  /. 

(3.44) 

Because 


plih)  =  -^-  and 
Piih)  =  + 

the  necessar>  condition  stated  in  equation  (3.5)  requires  that  for  times  /  =  /,  and  t  =  ry 
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Figure  15.  Trajectories  for  Uc  =  0;  From  (Ref.  8  :  p.  280] 


0  =  1  +  ;.  I  «a(/)  I  +  r* (0-^2(0  +  PalOWaCO  (3-45) 

Substituting  equations  (3.28)  and  (3.43)  into  equation  (3.45)  for  time  r,  yields 

0  =  1  +  /  I  u,(/,)  I  +  c,.r;(/,)  +  (  -  ;.)«a(/i)  (3.46) 

Substituting  equations  (3.28)  and  (3.44)  into  equation  (3.45)  for  time  tj  yields 

0  =  1  +  /  I  Ugiti)  I  +  ‘^l-''2(^2)  +  (  +  (3-47) 

Substituting  u{t)  =  0  into  equations  (3.46)  for  i  =  ii  and  equation  (3.47)  for  t  =  tj  yields 
-I=q^2'(fi)  (3.48) 

-l=c,.r;(/,)  (3.49) 

Now  solve  equations  (3.40)  tlirough  (3.49)  for  4r|(/,)  .  Equations  (3.48)  and  (3.49)  imply 
-Y2*('i)  =  -^2*(/2)  (3-50) 

Subtracting  equation  (3.44)  from  equation  (3.43)  results  in 


Figure  16.  Combined  Trajectories  for  U(t)=0  and  U(t)  =  (  +  -)Uc;  From  [Ref.  8; 

p.  281] 


Substituting  for  c,  from  equation  (3. 48)  into  equation  (3.51) 

C^2  ~  =  2/.V2(t|) 

Substitute  [/j  —  t,]  from  equation  (3.52)  into  equation  (3.42)  to  yield 


Equation  (3.40)  also  defines  .vK/j).  so  equating  it  with  equation  (3.53)  yields 


=  +  2/.v2*^(fi) 


Substituting  for  from  equation  (3.51)  into  equation  (3.54) 


(3.51) 

(3.52) 

(3.53) 

(3.54) 
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ar.'Vi) 


T777-7  =  -^i(^i)  +  2;.x2  (/,) 


Combining  terms  and  solving  for  jc,‘(ti)  in  equation  (3.55)  yields  ihe  equations  making  up 
the  optimal  control  law  for  a  weighted-time-fuel  control  response 


.  l+4/«^(/i)  ,2/.  s 


(3.56) 


»  I  +  4/U,(/i  ’) 


(3.57) 


and  equation  (3.40)  for  minimum-time  control  completes  the  optimal  control  law  set  of 
equations. 

The  trajectories  corresponding  to  the  optimal  control  law  define  the  switching  curves 
for  weight  _u-time-fuel  optimal  performance  and  are  shown  for  various  /.  values  and  ini¬ 
tial  conditions  in  Figure  17  on  page  37.  The  weighted-time-fuel  optimal  trajectories 
shown  in  Figure  17  on  page  37  also  illustrate  the  bang-dead  zone-bang  characteristic  of 
this  type  of  control. 

Using  the  switching  function  again,  equations  (3.56)  and  (3.57)  can  be  r:presented 
by 


5(.v(t))  S.Vi(t,)  - 


1  -f  4/i/^f/i) 


-V2(ri)  I  .^2(^1)  1 


(3.5S) 


The  switching  functions  derived  in  this  chapter  will  be  used  to  implement  the  optimal 
control  law  into  the  simulation  program  discussed  in  the  next  chapter. 

C.  DEADBAND  CONTROL  ABOUT  THE  ORIGIN 

A  control  deadband  about  the  origin  is  needed  in  addition  to  the  weighted-time-fuel 
optimal  control  law  and  corresponding  switching  curves  to  implement  attitude  control 
for  the  CER.  The  deadband  allows  for  the  CER  attitude  to  be  maintained  within  a  limit 
cycle  bounded  by  a  displacement  of  -I-/  —  1.25’  and  a  a  rate  of  +1  —  0.05’ In  addition 
to  the  limit  cycle  deadband,  an  inner  deadband  bounded  by  a  displacement  of 
-I-/  —  O.OF  and  a  rate  of  -H/  —  0.02’’"  is  installed  around  the  origin.  When  the  attitude 
controller  drives  the  system  to  within  this  inner,  or  zero  deadband,  the  simulation  or 
actual  logic  software  shuts  off  the  control  acceleration.  Without  deadbands  about  the 
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Figure  17,  Suitcliing  Curves  for  Weighted-Time- Fuel  Optimal  Performance;  After 
[Ref.  8  ;  p.  283] 

origin,  the  controller  would  chatter  back  and  forth  around  the  origin,  cycling  between 
positive  and  negative  control  accelerations.  The  residual  velocity  or  any  external  dis¬ 
turbance  torques  will  eventually  drive  the  system  back  out  to  the  outer  or  limit  cycle 
deadband  and  the  control  process  repeats.  An  illustration  of  the  two-deadband  control 
scheme  is  shown  in  Figure  18  on  page  38. 
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Figure  18.  Deadband  Control  around  the  Origin. 

The  simutcuion  step-size  has  a  major  elTect  on  setting  tlie  boundar}’  limits  for  the  in¬ 
ner  or  zero  deadband  around  the  origin.  If  the  step-size  is  larger  than  the  boundaries, 
it  completely  steps  over  the  zero  deadband  and  chatters  back  and  forth,  cycling  across 
the  deadband  without  shutting  off.  For  actual  operation  of  the  CER,  the  simulation 
step-size  could  be  thought  of  as  the  minimum  impulse  specilication  of  the  jet  thruster. 
If  the  minimum  impulse  time  (corresponding  to  the  step-size)  is  greater  than  the  zero 
or  inner  deadband  limits,  then  the  controller  will  not  be  able  to  shut  down  the  thruster 
before  the  system  is  already  back  out  of  the  dead  zone.  The  rate  boundary’  limits  are 
more  critical  due  to  the  switching  curves  used.  The  minimum-time  curves  approach  the 
origin  parabolically  and  the  rate  changes  faster  than  the  displacement  as  the  system  is 
driven  along  the  switching  curve. 

D.  DISCUSSION  OF  CONTROL  ACCELERATIONS. 

The  normalized  control  accelerations  to  be  used  for  each  axis  are  generated  by  the 
applied  torques  for  each  axis  acting  on  the  principal  moments  of  inertia  for  each  axis 
and  are  listed  as  follows: 
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"ci')  =  =  C  /„]■'  r,  -  f  r,  (3.59) 

u,y{t)  =  cb,.  =  C4^]-'  7;  =  y  (3.60) 

=  C4r’  7-.  =  j-  r,  (3.61) 


The  normalized  control  accelerations  are  also  components  of  the  total  amount  of 
acceleration  acting  on  all  three  axes,  which  includes  the  cross-coupling  effects  between 
axes  due  to  the  cross-products  of  inertia.  The  sum  of  the  accelerations  acting  on  each 
of  the  axes  are  listed  in  equations  (2.20),  (2.21),  and  (2.22),  and  are  listed  below  with  the 
control  accelerations  displayed  in  bold  face: 

((j  =  )^,.AT,-Ar,  +  ^r,  (3.62) 

+  (3.63) 

,0,s  =  ).t3  =  4r,y  r,  + jr.  (3.6J) 


The  bold  face  control  accelerations  above  are  the  only  components  used  in  the  optimal 
control  laws  derived  in  this  chapter.  Using  only  the  acceleration  components  generated 
by  the  principal  moments  of  inertia  to  drive  the  respective  states  (angular  position  and 
velocity),  of  each  axis  to  zero,  decouples  each  axis  attitude  control  from  the  other  axes. 
It  will  be  shown  in  the  next  chapter  that,  if  the  sum  of  the  acceleration  components  (for 
each  axis)  generated  by  the  cross-products  of  inertia  is  greater  than  the  control  acceler¬ 
ation  component  for  that  axis,  then  the  final  state  of  zero  is  never  reached. 
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IV.  SIMULATION  AND  RESULTS 


A.  SIMULATION  PROGRAM 

Attitude  control  performance  for  the  CER  baseline  configuration  is  computed 
through  computer  simulation  of  a  model  describing  the  CER  in  terms  of  its  equations 
of  motion,  total  moment  of  inertia  tensor,  summation  of  applied  torques,  selected  initial 
conditions,  and  optimal-fuel  weighting  factors.  The  equations  of  motion  used  in  the 
simulation  are  represented  by  equation  (2.17),  which  is  presented  in  component  form  by 
equations  (3.62),  (3.63),  and  (3.64).  The  inverse  moment  of  inertia  tensor  values  in¬ 
cluded  in  these  equations  were  calculated  by  substituting  the  moment  of  inertia  compo¬ 
nents  for  each  target  capture  case  listed  in  Figure  9  on  page  19  into  equation  (2.19).  The 
baseline  configuration  control  torques  included  in  the  equations  are  listed  in  Section  1.5 
on  page  4. 

1.  Simulation  Block  Model 

The  simulation  program  was  written  using  the  TUTSI M  computer  simulation 
language  to  accomplish  a  dynamic  simulation  of  the  block,  model  shown  m  Figure  19 
on  page  41  (Ref  9;  pp.  1. 1-1.3).  Not''  that  the  angular  acceleration  components  are 
integrated  twice  in  succession  to  yield  angular  velocity  and  angular  position  for  each 
axis.  The  switching  curves  and  the  deadbands  about  the  origin  described  in  Chapter 
three  are  used  to  determine  the  control  acceleration  required  to  drive  the  system  to  the 
desired  final  state  of  zero  for  each  axis.  The  angular  velocity  and  position  yielded  by  the 
integrators  define  the  location  of  the  system  on  the  state-space  plot  for  each  axis  as 
previously  shown  in  Figure  17  on  page  37.  The  control  acceleration  components  ap¬ 
plied  to  drive  the  system  to  zero  help  define  the  switching  curves  and  are  listed  bold  face 
in  equations  (3.62),  (3.63),  and  (3.64).  The  acceleration  components  (for  each  axis) 
generated  by  the  cross-products  of  inertia  act  as  disturbance  accelerations  driving  the 
system.  They  are  input  through  the  integrators  and  subsequently  affect  the  systems  lo¬ 
cation  on  the  state  space  plots.  From  Figure  19  on  page  41  it  is  easy  to  see  that  if  the 
sum  of  the  acceleration  components  for  an  axis  generated  by  the  cross-products  of  in¬ 
ertia  is  greater  than  the  control  acceleration  for  that  axis,  then  the  final  state  of  zero 
may  not  be  reached.  The  simulation  runs  displaying  unstable  control  are  expected  to 
have  control  accelerations  less  than  the  cross-product  generated  acceleration  term  sums. 


Figure  19.  Simulation  Clock  Model  of  the  CER 


2.  Simulation  Inputs 

Ihe  simulation  program  is  listed  in  “AI’PEXDIX  A.  TUISI.VI  SIMU¬ 
LATION  PROGIUAM"  on  page  .58  and  requires  the  following  inputs  to  dynamically 
simulate  three-axis  attitude  control  for  the  CER: 

•  riirec  control  torques  (fi-lhfy. 
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(Block  73) 

Yr,-^  (Block  249) 

Yr.- (Block  450  > 

•  Nine  inverse  moment  of  inertia  values  C/D'* 

y->  (Block  70) 

y->  (Block  240) 

(Block  439) 

J 

(Block  6S) 

--(Block  239) 
f 

y- (Block  443) 

y- (Block  71) 

—  -(Block  243) 

./ 

y- (Block  44) 

•  Three  lambda  values  (/  =  0  for  minimum-time): 

—  (Block  15) 

/ ^  —  (Block  213) 

/ .  — ►  ( Block  41.'') 

•  Initial  conditions  for  each  axis  of  angular  position  in  radians  and  angular  velocity 
in  radians  second; 

Roll  position  —  (Block  1) 

Roll  velocity  —  (Block  2) 

Pitch  position  —  (Block  201) 

Pitch  velocity  —  (Block  200) 

'^'aw  position  —  (Block  401) 

■^'aw  velocity  —  (Block  400) 

A  detailed  description  of  the  different  types  of  blocks  used  in  the  simulation  program  can 
be  found  in  Ref  9;  pp.  6.1-6.83. 

3.  Simulation  Cases 

The  CER  baseline  configuration  was  simulated  without  a  target  and  with  a  tar¬ 
get  captured  in  difTerent  net  locations  as  shown  previously  in  Figure  9  on  page  19.  The 
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radian  equivalents  for  the  initial  conditions  of  2’  angular  position  and  0.2’  sec  angular 
velocity  are  assigned  to  each  axis  for  all  the  simulation  runs  for  comparison  to  allow  for 
easy  \  isualization  of  the  dynamic  movement.  The  initial  conditions  correspond  to  the 
system  state  values  after  being  oO'set  from  the  desired  reference  frame  by  some  disturb¬ 
ance  impulse.  The  angular  velocities  were  assumed  to  be  small  and  the  angular  velocity 
products  were  dropped  from  equations  (2.13),  (2.14),  and  (2.15).  The  simulation  model 
is  now  only  applicable  to  manuevers  using  small  angular  velocities.  The  model  is  ideal 
for  autonomous  attitude-hold  simulation  since  only  small  position  and  velocity  initial 
condition  offsets  arc  used.  The  model  is  not  applicable  to  simulation  of  commanded 
rotations  (slewing)  for  minimum-time  control  of  the  CER  other  than  making  small  angle 
position  adjustments  since  the  corresponding  maximum  angular  velocities  attained  dur¬ 
ing  the  manuevers  are  too  large.  A  weighted-time-fuel  optimal  control  could  be  used  for 
slewing  ntanuevers  since  the  maximum  angular  velocities  attained  would  be  minimized 
by  the  dcadzonc  characteristic  of  the  control.  This  results  in  a  slower  time  response  to 
drive  the  states  to  zero.  A  more  detailed  discussion  will  follow  later  in  this  chapter. 

To  determine  control  stability  for  the  selected  target  capture  cases,  only 
minimum-time  simulations  were  initially  run.  The  optimal-fuel  weighting  fitetor  lambda 
(/. )  for  each  axis  was  according!}'  assigned  a  zero  value. 

B.  B.\SELI.NE  CONFIGURATION  SIMULATION  RESULTS 

The  angular  position  and  velocity  for  each  a.xis  were  driven  to  zero  from  the  selected 
initial  conditions  using  minimum-time  control.  Of  main  interest  is  the  maximum  over- 
slh'Oi  attained  by  the  states  while  being  driven  to  zero  and  the  time  it  took  to  initially 
arrive  near  zero. 

Th.e  CER  baseline  configuration  using  1.0  Ibf  thrusters  was  simulated  for  all  target- 
capture  cases  and  the  results  are  shown  in  Figure  20  on  page  44  along  with  the  moment 
of  inertia  tensors  corresponding  to  each  case.  The  left-handed  arrows  indicate  the  axes 
where  the  desired  final  state  of  zero  was  never  reached  and  the  unstable  overshoots  are 
indicated  by  the  infinity  (oo)  symbols.  The  five  case  a.xis  combinations  experiencing 
unstable  control  are  from  the  850  pound  point  mass  target  capture  cases.  The  absolute 
value  of  the  control  acceleration  component  for  each  axis  experiencing  unstable  control 
is  less  than  the  absolute  value  of  the  sum  of  the  acceleration  components  generated  from 
the  cro'^s-products  of  inertia. 
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CASE 

MAXOVERSH<X)T 

FROM  I.C.  =  2°/0.2°^”' 

&  TIME  TO  ZERO  (sec) 

NO  TARGET 

R 

2.1  /  2.9  /  1.5 

P 

2.1  /  2.3  /  1.5 

Y 

2.1  /  2.9  /  1.5 

CASE1A  (Xmax) 

R 

2.1  /  2.5  /  1.4 

P 

2.2  /  0.6  /  8.3 

Y 

oo  ^ 

CASEIB(Ymax) 

R 

2.0  /  2.1  /  1.9 

P 

OO  ^ - 

Y 

OO  _ 

CASE1C(2max) 

R 

2.5  /  0.5  /  15.6 

P 

2.2  /  0.7  /  6.6 

Y 

2.1  /  1.2  /  3.8 

CASE1D 

R 

2.0  /  2.0  /  2.2 

(X=Z.Y=^1) 

P 

OO  S  ■  ■ 

Y 

oo  / 

CASE2 

R 

2.0  /  1.9  /  2.2 

(MAN+MMU) 

P 

2.1  /  0.9  /  4.9 

Y 

2.1  /  0.8  /  6.6 

MOMENT  OF  INERTIA 
TENSOR 


Figure  2(1.  Baseline  Configuration  Simulation  Results 


1  he  "No  'I'argct"  case  and  "Case  lA"  angular  acceleration  components  were  hand- 
calculated  using  equations  (3.62),  (3.63).  and  (3.64)  to  gain  further  insight  into  the  sim¬ 
ulation.  fhe  control  acceleration  magnitude  is  shown  bold  face  as  follows: 

•  No  1  argot: 
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Similar  calculations  for  the  rest  of  the  cases  were  completed  by  the  author  and  for  every 
case  of  unstable  control,  the  absolute  value  of  the  control  acceleration  magnitude  was 
less  than  the  absolute  value  of  the  sum  of  the  acceleration  components  generated  from 
the  cross-products  of  inertia. 

The  maximum  overshoot  values  displayed  in  Figure  20  on  page  44  were  taken  from 
the  simulation  runs  in  state  space  plot  form  for  each  of  the  axes  of  every  case.  The  times 
to  initially  arrive  to  zero  were  taken  ofT  the  simulation  runs  of  the  angular  positions 
ver'us  time  plots  for  every  case.  The  simulation  plots  are  found  in  “APPENDIX 
B.  SIML'L.ATION  PLOTS"  on  page  62  and  the  angular  positions  and  velocities  are 
displayed  in  radians  and  radians  second  respectively.  Note  that  some  of  the  simulations 
were  allowed  to  run  long  enough  to  show  the  limit  cycle  resulting  from  the  inner  and 
outer  deadband  regions. 
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The  author  tried  various  combinations  of  torque  summation  values  for  the  three 
axes  in  an  attempt  to  reduce  the  number  of  unstable  cases  without  changing  the  CER 
structure  or  thruster  location.  The  best  combination  found  after  many  iterations  was 


^r,  =  34  (ft-ibo 

(4.1) 

^7^,=  100  (ft-lbO 

(4.2) 

100  (ft-lbO 

(4.3) 

which  yielded  unstable  control  about  four  case/axis  combinations  as  shown  in 
Figure  21  on  page  47.  The  simulation  plots  are  not  displayed  in  Appendix  B  to  avoid 
confusion. 

C.  PROPOSED  CONFIGURATION  AND  SIMULATION  RESULTS 
1.  Proposed  Configuration 

The  cross-products  of  inertia  must  be  reduced  and  the  summation  of  control 
torques  about  the  axes  may  have  to  be  increased  to  achieve  stable  control  about  all  axes 
of  every  target  capture  case.  The  target  capture  mechanism  can  be  raised  up  to  center 
on  the  X  axis  to  reduce  the  cross-products  of  inertia.  1  his  enables  the  target  to  be  cap¬ 
tured  closer  to  the  CER  center  of  gravity,  resulting  in  a  reduction  of  the  magnitude  of 
the  vector  fj  and  the  corresponding  skew-synunetric  matrix  as  shown  previously  in  Fig¬ 
ure  8  on  page  18.  Equation  (2.36)  shows  that  a  reduction  in  the  skew-synunetric  matrix 
reduces  the  total  moment  of  inertia  tensor.  The  following  thruster  sizes  are  proposed  to 
be  used  at  the  original  baseline  locations  since  the  torque  combination  from  equations 
(4.1),  (4.2),  and  (4.3)  reduced  the  number  of  unstable  control  axes  for  the  baseline  con¬ 
figuration: 


=  11.3  Ibf 

(4.4) 

F^  =  33.33  Ibf 

(4.5) 

F,  =  25  ibf 

(4.6) 
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CASE 

MAX  OVERSHOOT 

o  °/ssc 

FROM  I.C.  =  2  /0.2 
&TIME  TO  ZERO  (sec) 

NO  TARGET 

R 

2.0  /  9.9  /  0.3 

P 

2.0  /  14.3  /  0.3 

Y 

2.0  /  14.3  /  0.3 

CASE1A  (Xmax) 

R 

2.0  /  8.0  /  0.5 

P 

2.0  /  2.6  /  1.2 

Y 

2.1  /  2.0  /  2.1 

CASE1  B(Ymax) 

R 

oo 

P 

2.0  /  2.3  /  2.0 

Y 

2.0  /  2.3  /  2.0 

CASEIC(Zmax) 

R 

OO 

P 

2.0  /  3.9  /  1.0 

Y 

2.0  /  8.0  /  0.45 

CASE1D 

(X=Z.Y=1) 

R 

OO 

P 

oo  ^ - 

Y 

2.0  /  3.8  /  1.1 

CASE2 

R 

2.0  /  4.8  /  0.8 

(MAN+MMU) 

P 

2.0  /  4.7  /  0.8 

Y 

2.0  /  4.9  /  0.8 

MOMENT  OF  INERTIA 
TENSOR 


Figure  21.  Baseline  Configuration  Simulation  Results  With  T.\  =  34,  Ty  =  Tz=  100 
(ft-lbO 


Ihc  CER  proposed  configuration  is  illustrated  in  Figure  22  on  page  48,  and 
shows  an  example  of  a  thruster  pair  liring  to  produce  a  positive  roll  about  A'.  Note  the 
capture  mechanism  is  now  attached  so  as  to  be  centered  on  A' . 
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2.  Proposed  O  ’'figuration  Analysis  and  Simulation  Results 

The  methods  described  in  Chapter  Two  for  computing  the  moment  of  inertia 
tensors  for  the  selected  target  capture  cases  were  used  to  generate  the  tensors  for  the 
proposed  configuration.  The  CER  proposed  configuration  was  simulated  for  all  target 
capture  cases  and  the  results  arc  shown  in  Figure  23  on  page  49  along  with  the  moment 
of  inertia  tensors  corresponding  to  each  case.  The  left-handed  arrow  indicates  unstable 
control  about  the  X  a.xis  for  Case  1C  target  capture.  All  the  other  target  capture  cases 
have  stable  control  with  satisfactory  ma.ximuni  overshoot  values. 

Case  1C  is  the  850  pound  point  mass  captured  at  the  bottom  of  the  capture  net. 
The  target  was  brought  closer  to  tlic  center  of  the  net  until  the  simulation  produced 
stable  control.  Satisfactory  stable  control  occurred  with  the  target  located  +  2  feet  from 
the  net  center  in  the  Z  direction,  and  is  represented  by  Case  IC.l  in  Figure  23  on  page 
49. 


48 


Figure  23.  CER  Proposed  Configuration  Simulation  Results 


All  of  the  moment  of  inertia  tensors  in  Figure  23  on  page  49  have  smaller 
component  magnitudes  than  the  baseline  configuration  tensors  shown  in  Figure  20  on 
page  44.  I  he  time  response  to  initially  drive  the  system  to  zero  is  satisfactory.  I  he 
simulation  plots  for  the  proposed  configuration  are  shown  in  Appcndi.x  B. 
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3.  Proposed  Configuration  With  Baseline  I  hrusters 

The  proposed  configuration  was  simulated  using  1.0  Ibf  thrusters  in  the  origijial 
baseline  locations  to  minimize  changes  to  the  baseline  configuration.  The  simulation 
results  are  shown  in  Figure  24  on  page  51  and  the  left-handed  arrows  indicate  unstable 
control  about  2  axes  of  Case  ID  and  1  axis  of  Case  IB.  1  he  simulation  plots  are  not 
included  in  Appendix  Bto  avoid  confusion. 

D.  FUEL  OPTIMAL  ANALYSIS  FOR  THE  PROPOSED  CONFIGURATION 
llie  proposed  configuration  was  simulated  using  the  vveighted-time-fuel  optimal 
switching  curves  with  a  weighting  factor  ).  =  100  for  each  axis.  The  No  Target  case  was 
simulated  and  the  simulation  plots  can  be  seen  in  Appendix  B.  The  simulation  plots  for 

A 

a  positive  roll  about  the  X  axis  for  minimum-time  and  weighted-time-fuel  optimal  con¬ 
trol  are  compared  in  Figure  25  on  page  52.  The  deadzone  on  the  optimal  control  plot 
minimizes  the  maximum  angular  velocity  value  attained.  The  time  response  to  initially 
drive  the  system  to  zero  is  increased  when  compared  to  minimum-time  control. 

The  simulation  plots  of  position  versus  time  for  minimum-time  and  optimal  control 
with  ).  =  lUO  give  the  time  to  complete  one  limit  cycle  and  the  thruster  "on  times"  for 
one  limit  cycle  as  follows: 

•  Time  to  complete  one  limit  cycle  starting  from  origin  for  2=0  (nrinimum-time); 

Roll  ->  90  sec 
Pitch  -►  58  sec 
Yaw  -♦  58  sec 

•  Time  to  complete  one  limit  cycle  starting  from  origin  for  X  =  100  (optimal-fucfi: 

Roll  140  sec 
Pilch  -♦110  sec 
Yaw  -♦  1 10  sec 

•  "On  times"  for  one  limit  cycle  with  X  =  0: 

Roll  -♦  3  sec 
Pitch  -♦  2  sec 
Yaw  ->  2  sec 

•  "On  times"  for  one  limit  cycle  with  X  =  100: 

Roll  -♦  I  sec 
Pitch  -♦  0.4  sec 
Yaw  -♦  0.4  sec 

The  number  of  cycles  completed  during  30  minutes  (1800  seconds)  of  operation  is: 

•  For  X  =  0: 
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VlMHrUT  x2  INITIAL  CONDITION 
8.1745006  ^  ^ 

8 . 1390000 

8 . 1047000 

0.0090000 

0.0349000 


NO  TARGET  ROLL  PROPOSED 


-8.0349000 

-0.0090000 

-0.1047000 

-0.1390000 


-0.1745000 


-0.007200 INPUT  xl  INITIAL  CONDITION 


0.0072000 


VIMHMIT  x2  INITIAL  CONDITION 
0.0072000  ^ 

0.0090000 

0.dS235C.n 

0.0349040 

0.0174520 


-0.0174520 


NO  TRGT  ROLL  PROPOSED  OPT.  (L=1O0) 


-0.0349040 


-0.0523500 


-0.0090000 

-0.0072000 


-0.007200 INPUT  xl  INITIAL  CONDITION 


0.0072000 


Figure  25.  Comparison  of  Minimum-Tiine  and  Optimal  Fuel  (Lambda  =  100)  Roll 


AaIs  Control  for  "No  Target" 


r.  „  1800  sec  , 

Roll  -►  -r-; - —  =  12.9  cvcles 

1-lU  see; cycle 

Pitch  and  Yaw  -»  ^  _  j(3  4  cvcles 

1  lU  sec.'cvcle 


Total  "on  time"  for  one  thruster  is  found  by  multiplying  the  number  of  cycles  by  the  on 
time  per  cycle  to  yield 


•  For  /  =  0: 


Roll  -+  60  sec 

Pitch  and  ^'a\v  -♦  62  sec 

•  For  /  =  100; 

Roll  ->  12.9  sec 
Pitch  and  Yaw  -»  6.6  sec 

Two  thrusters  fire  per  axis  or  six  thrusters  total  fire  as  needed  for  three-axis  control. 
The  thruster  flow  rate  is  chosen  for  analysis  to  be  0.5  pounds,' second  of  nitrogen  (Aj). 
Fuel-optimal  analysis  for  the  proposed  configuration  is  summarized  in  Figure  26  on 
page  54. 

E.  END-OF-xMISSION  ANALYSIS  FOR  THE  PROPOSED  CONFIGURATION 

The  CER  mass  will  decrease  by  approximately  150  pounds  at  the  end  of  its  mission 
when  the  total  nitrogen  propellant  is  used  up.  The  CER  mass  is  reduced  to  700  pounds 
and  when  a  worst-case  850  pound  target  is  captured  the  center  of  gravity  is  closer  to  the 
target  than  at  the  beginning  of  the  mission.  The  shift  in  the  center  of  gravity  does  not 
effect  the  summation  of  torques  about  the  axes  since  the  thrusters  are  fired  in  pairs.  The 
decrease  in  the  CER  mass  (.V/,)  reduces  the  total  moment  of  inertia  tensors  for  all  target 
capture  cases  as  shown  in  equation  (2.36).  The  moment  of  inertia  tensor  for  the  CER 
without  a  target  at  the  end  of  the  mission  is  reduced  to: 


^CER  — 


32.6  0  0 

0  45  (J 

0  0  45 


and  the  reduction  improves  minimum-time  control.  The  reduction  of  the  moment  of 
inertia  tensors  improves  the  minimum-time  control  for  all  target  capture  cases  as  shown 
in  Figure  27  on  page  55. 
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FUEL  OPTIMAL  ANALYSIS 
FOR  PROPOSED  CONFIG. 

OF  CER  WITHOUT  TARGET 

LAMBDA  =  0 
(MIN  TIME) 

LAMBDA  =  100 
(OPT.  FUEL) 

TOTAL  "ON"  TIME  DURING 
30  MINUTES  OPERATION 
WITH  6  THRUSTERS 

368  SECONDS 

52 SECONDS 

TOTAL  N2  USED  DURING 

30  MINUTES  OPERATION 

18.4  LBS 

2.6  LBS 

"TIME  TO  ZERO"  FROM  I.C.: 

O  o/  SEC 

2  !  Q.2  ' 

Bl 

p  J  3  SECONDS 

Y  ^ 

Figure  26.  Fuel  Optimal  Analysis  for  Proposed  ConFiguration 
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Figure  27.  End-of-Mission  Siimilntioii  Results  lor  the  Proposed  Configuration 


V.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

The  CER  baseline  configuration  of  1.0  Ibf  thrusters  and  capture  mechanism  at¬ 
tached  at  the  bottom  of  the  CER  did  not  provide  effective  control  during  capture  of  a 
worst-case  850  pound  target.  Changing  the  thruster  size  and  location  to  increase  the 
sunimation  of  torques  about  the  axes  did  not  provide  effective  control  during  target 
capture.  The  unstable  control  cases  were  shown  to  be  caused  by  the  acceleration  com¬ 
ponents  generated  from  the  cross-products  of  inertia. 

The  CER  proposed  configuration  with  increased  summation  of  torques  about  the 
axes  and  capture  mechanism  attached  along  the  X  axis  provides  an  effective  control 
during  all  cases  of  target  capture  e.xcept  for  Case  1C  (  +  -  roll).  The  unstable  control 
became  stable  when  the  target  distance  from  the  center  of  the  net  was  reduced  from  3.5 
feet  to  2  feet  on  the  Z  axis.  Control  stability  was  dramatically  improved  by  moving  the 
capture  mechanism  closer  to  the  CER  center  of  gravity  resulting  in  a  reduction  of  the 
cross-products  of  inertia  for  target  capture  operation. 

Weighted-time-fuel  optimal  control  with  /  =  100  reduced  propellant  use  by  S5%  for 
the  selected  example  but  increased  the  time  response  to  drive  the  system  to  zero  by  a 
factor  of  li).  The  maximum  angular  velocity  values  attained  are  minimized  by  the  dead 
zone  of  the  optimal  control.  Weighted-time-fuel  optimal  control  could  be  used  during 
coasting  periods  between  orbits  when  a  minimum-time  response  is  not  required  or  for 
slow  commanded  rotations  (slewing)  since  the  angular  velocities  are  minimized. 

Control  performance  is  improved  at  the  end  of  the  mission  due  to  the  decrease  in 
mass  of  the  CER.  The  center  of  gravity  shift,  as  the  mass  of  the  CER  is  reduced,  does 
not  efiect  the  summation  of  control  torques  about  the  axes  for  the  CER  with  target. 

B.  RECO.\TMENDATIONS 

The  capture  mechanism  should  be  attached  along  the  X  axis  to  reduce  the  distance 
between  the  CER  center  of  gravity  and  the  target  center  of  gravity.  The  cross-products 
of  inertia  will  be  reduced  and  control  stability  will  be  improsed. 

The  baseline  thruster  configuration  should  be  changed  to  provide  control  torques 
as  follows  if  the  CER  model  used  is  similar  to  the  model  presented  in  this  research: 


y  r,=  34  ft-lbf 
V  T,  =  100  ft-lbf 

The  baseline  thruster  location  could  be  used  to  facilitate  storage  of  the  CER  and  to 
minimize  changes  in  the  baseline  configuration  if  the  thruster  size  is  increased  to 

=  11.3  Ibf 
F,  =  33.3  Ibf 
R  =  25  Ibf 

Weighted-time-fuel  optimal  control  could  be  used  during  coasting  periods,  between 
orbits  or  for  slow  or  small  commanded  rotations,  when  a  minimum-time  response  is  not 
required.  Operation  of  the  CER  with  target  may  require  immediate  minimum-time  re¬ 
sponse  and  probably  should  not  use  weighted-time-fuel  optimal  control. 

The  proposed  attitude  control  scheme  is  simple  and  elTective  but  is  sensitive  to  tar¬ 
get  capture  location  in  the  capture  net.  Ensuring  the  target  is  captured  as  close  as  pos¬ 
sible  to  the  CER  center  of  gravity  will  reduce  the  cross-products  of  inertia  and  improve 
the  design  and  operational  attitude  control  performance. 

The  control  laws  used  to  implement  simulation  of  the  CER  model  are  in  terms  of 
the  control  accelerations  generated  from  the  principal  moments  of  inertia  and  of  the 
angular  position  and  velocity  for  each  of  the  axes.  A  more  complicated  control  law  could 
be  used  in  terms  of  the  control  accelerations,  angular  positions,  and  angular  velocities 
for  the  three  axes  in  combination.  ,A  switching  surface  or  volume  corresponding  to  the 
control  law  would  be  less  sensitive  to  the  cross-product  of  inertia  elTccts  as  opposed  to 
the  simple  switching  curves  used  in  these  analyses. 

The  research  presented  should  assist  the  CER  preliminarx’  design  phase  team  in 
completing  some  of  the  attitude  control  trade-off  analyses.  Ultimately,  this  will  result 
in  a  more  practical  and  safe  crew  and  equipment  retriever  for  space  station. 
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APPENDIX  A.  TUTSIM  SIMULATION  PROGRAM 


The  TUTSIM  simulation  program  will  simulate  three-axis  stabilization  for  a  model 
described  in  terms  of  the  sum  of  the  torques  applied  and  the  system  moment  of  inertia 
tensor  including  non-zero  cross-product  terms.  The  program  is  shown  in  Figure  28  on 
page  59,  and  continued  in  Figure  29  on  page  60.  and  Figure  30  on  page  61.  The  first 
199  block  numbers  are  reserved  for  the  Roll  axis,  blocks  200  through  399  are  for  the 
Pitch  axis,  and  blocks  400  through  599  are  for  the  Yaw  a.xis.  The  interconnecting  inverse 
moment  of  inertia  terms  are  interspersed  throughout  the  program.  The  inputs  to  the 
simulation  program,  are  listed  in  Section  IV. 2  on  page  41.  The  simulation  is  sensitive  to 
the  step-size  used  since  the  inner  or  zero  deadband  boundrics  are  small.  A  step-size  of 
.001  worked  nicely  most  of  the  time.  If  the  simulation  chatters  at  the  origin,  reduce  the 
step-size.  The  specific  example  to  follow  is  for  the  CUR  baseline  configuration  with  no 
target. 
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I'ROFESSIONAL  VERSION  OF  TUTSIH 


Model  File:  cerbsnt 

Dale:  2  /  28  /  1989 

Time:  13  :  45 

Timing:  0.0400000  .DELTA  ;  lO.OOOEiOT  ,Ri.«CE 

PlotBlocks  and  Scales: 

Fo  rma  t ; 


BlockNo , 

Plot-MINlmum, 

Plot-MAXlraura 

Comment 

Ho  17 

1  , 

-  5 . 0000 

5 . 0000 

INPUT  xl 

INITIAL 

CONDITION 

VI 

Y2 

Y3 

Y', 

2  . 

-5.0000 

5 . 0000 

INPUT  x2 

INITIAL 

CONDITION 

0.0350000 

1  INT 

2 

; INPUT  xl  INITIAL  CONDITION 

0.0035000 

2  INT 

67 

INPUT  x2  INITIAL  CONDITION 

3  SUM 

1 

10 

xl+(0.5/Ucx)x21x2) 

0.0000 

4  REL 

8 

-8 

-8 

slgn(xl) 

0.5000000 

5  CON 

1 .0000 

6  CON 

0.0000 

7  CON 

8  HUL 

70 

73 

(n/j)Tx-Ucx 

to  HUL 

2 

11 

66  ; 

(0.5/Ucx)x21x2| 

11  ARS 

2 

1x2) 

12  SUM 

1 

13 

xl»((l+4LUcx)/(2Ucx))x21x2) 

13  HUL 

2 

11 

22  ; 

((U4LUcx)/(2Ucx))x21x21 

4  0000 

14  CON 

0.0000 

15  CON 

INPUT  lambda ."L) , (optimal  wtng) 

16  HUL 

14 

15 

8  : 

4LUcx 

17  SUM 

6 

16 

li4LUcx 

2.0000 

18  CON 

19  HUL 

18 

8 

2Ucx 

20  IFE 

44 

47 

21 

begin  dead  zone 

21  IFE 

43 

7 

-47 

leave  dead  zone 

22  DIV 

17 

19 

(l+4LUcx)/(2Ucx) 

37  INV 

38 

true  If  outside  zero  rect. 

38  AND 

48 

41 

Inside  zero  rect. 

350.000E-06 

39  CON 

INPUT  zero  rect.  velocity  (x2) 

40  ARS 

1 

1x1) 

41  SUM 

39 

-11 

0.001- (x2J 

42  IFE 

38 

7 

46 

inside  zero  rect  .  7 

4  3  HUL 

4 

3 

(xl+(0. 5/Ucx)x2 (x2 ] )slgn(xl) 

44  HUL 

4 

12 

(Block  «12)slgn(xl) 

45  HUL 

1 

2 

xl*x2 

46  IFE 

45 

.4 

20 

set  U  for  Ist6i3rd  quad 

0 . 0000 

4  7  REL 

8 

O 

-8 

-8 

slgn(x2) 

48  SUM 

c 

49 

-40 

0.01- 1x1 1 

i;5.000E-06 

49  CON 

INPUT  zero  rect.  position  (xl) 

56  ORR 

59 

61 

outside  outer  rect.  7 

57  INV 

56 

true  If  Inside  outer  rect. 

0.0218000 

58  CON 

INPUT  outer  rect.  position 

59  SUM 

-58 

40 

1x1 1-1.25 

872 .000E06 

60  CON 

INPUT  outer  rect.  velocity 

61  SUM 

-60 

11 

1x2) - .01 
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62  IFE 

38 

7 

92 

U-0  If  Inside  zero  rect. 

1 . 0000 

63  SRS 

37 

57 

l.D.  and  loop  thru  regions 

6'i  AND 

63 

37 

roinpnre  loop  niul  locotlon 

65  IFE 

69 

62 

7 

U-0  until  outside 

66  DIV 

5 

8 

0. 5/Ucx 

67  SUM 

65 

291 

990 

sum  of  control  Inputs 

0 . 0000 

68  CON 

INPUT  d/j 

69  HUE 

-68 

73 

(-d/J)Tx 

0.02526/.0 

70  CON 

INPUT  a/J 

0.0000 

71  CON 

INPUT  g/J 

72  HUE 

71 

73 

(E/J)Tx 

3 .0000 

73  CON 

INPUT  Tx 

0.0035000 

200  INT 

292 

INPUT  y2  INITIAE  CONDITION 

0.0350000 

201  INT 

200 

INPUT  yl  INITIAE  CONDITION 

202  SUM 

201 

207 

yl+(0.5A)cy)y2iy2j 

0.0000 

203  REE 

209 

-209 

-209 

sign(yl) 

201 

209  MUE 

239 

299 

(e/J )Ty-Ucy 

0.5000000 

205  CON 

206  DIV 

205 

209 

0 . 5/Ucy 

207  MUE 

200 

208 

206 

(0.5/Ucy)y2|y2j 

208  ABS 

200 

[y21 

0.0000 

209  CON 

0.0000 

210  REE 

209 

-209 

-209 

slgn(y2) 

200 

211  IFE 

212 

209 

-210 

leave  deadzone 

212  HUE 

203 

202 

(yl+(0.5/Ucy)y2[y2)).slgn(yl) 

0.0000 

213  CON 

INPUT  lambcla'-l^,  (opt .  fuel  wtng) 

219  HUE 

200 

208 

298 

((H9EUcy)/(2Ucy))y2(y21 

215  SUM 

201 

219 

yl+((l+9EUcy)/(2Lcy))y2(y2] 

216  HUE 

203 

215 

(Block  «215)slgn(yl) 

217  IFE 

216 

210 

211 

begin  deadzone 

218  MUE 

201 

200 

yl*y2 

219  IFE 

218 

-203 

217 

set  Uy  for  lst63rd  quadrants 

220  IFE 

222 

209 

219 

Inside  zero  rectangle? 

221  IFE 

222 

209 

220 

Uy-0  If  Inside  zero  rectangle 

222  AND 

225 

223 

Inside  zero  rectangle? 

223  SUM 

229 

-208 

0.001-(v2) 

350.000E-0b 

229  CON 

INPUT  zero  rect.  velocity  (y2) 

225  SUM 

226 

-227 

0.01-(yl] 

175.000E-06 

226  CON 

INPUT  zero  rect.  position  (yl) 

227  ABS 

201 

lyil 

0.0218000 

228  CON 

INPUT  outer  rcctiingle  position 

872 .OOOE-06 

229  CON 

INPUT  outer  rectangle  velocity 

230  SUM 

-228 

227 

lyll-1.25 

231  SUM 

-229 

208 

ly2)-.oi 

232  ORR 

230 

231 

outside  outer  rectangle? 

233  INV 

232 

true  if  inside  outer  rectangle 

239  INV 

222 

true  if  out.slde  zero  rectangle 

1.0000 

235  SRS 

239 

233 

Identify  and  loop  thru  regions 

236  AND 

235 

239 

compare  loop  and  location 

237  IFE 

236 

221 

209 

Uy-O  until  outside 

0. 01818^3 

239  CON 

INPUT  (e/J) 

0.0000 

290  CON 

INPUT  (b/J) 

291  MUE 

-290 

299 

(-b/J)Ty 

292  SUM 

69 

237 

hUU 

sum  of  coritrol  torques 

0  0000 

293  CON 

INPUT  (h/J) 

299  MUE 

-293 

299 

(-h/J)Ty 

295  MUE 

19 

213 

209 

^LUcy 

296  SUM 

6 

295 

l+4LUcy 
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247  MUL 

18 

204 

2Ucy 

248  DIV 

246 

247 

(l44LUcy)/(2Ucy) 

3.0000 

249  CON 

INPUT  Ty 

0.0035000 

400  INT 

445 

INPUT  z2  INITIAL  CONDITION 

0.0350000 

401  INT 

400 

INPUT  zl  INITIAL  CONDITION 

402  SUM 

401 

407 

zl+(0.5/Ucz)z2|z2) 

0.0000 

403  REL 

404 

-404 

-404 

9lgn(zl) 

401 

404  MUL 

441 

450 

(l/j)Tz-Ucz 

0 . 5000000 

405  CON 

406  DIV 

405 

404 

0.5/Ucz 

407  HUL 

400 

408 

406 

(0.5/Ucz)z2lz21 

408  ABS 

400 

lz2) 

0.0000 

409  CON 

0 . 0000 

410  REL 

404 

-404 

-404 

slgn(z2) 

400 

411  IFE 

412 

409 

-410 

leave  deadzone 

412  MUL 

403 

402 

(z  1-4(0 . 5/Ucz)z21z2]  )slgn(zl) 

0 . 0000 

413  CON 

INPUT  lambda~L, (opt .  fuel  wtng) 

414  MUL 

400 

408 

449 

((l44LUcz)/(2Ucz))z2(z2) 

415  SUM 

401 

414 

zl4((l44LUcz)/(2Ucz))z2iz2J 

416  MUL 

403 

415 

(Block  W415)slgn(zl) 

417  IFE 

416 

410 

411 

begin  deadzone 

418  HUL 

401 

400 

zl*z2 

419  IFE 

418 

-403 

417 

set  Uz  for  lst&3rd  quadrants 

420  IFE 

422 

409 

419 

Inside  zero  rectangle? 

421  IFE 

422 

409 

420 

Uz-0  If  Inside  zero  rectangle 

422  AND 

425 

423 

inside  zero  rectangle? 

423  SUM 

424 

-408 

0.001-|z2) 

350.000E-06 

424  CON 

INPUT  zero  rect.  velocity  (z2) 

425  SUM 

426 

-427 

O.Ol-lzl) 

175.000E-06 

426  CON 

INPUT  zero  rect.  position  (zl) 

427  ABS 

401 

0.0218000 

428  CON 

INPUT  outer  rectangle  position 

872.000E-06 

429  CON 

INPUT  outer  rect.  velocity 

430  SUM 

•428 

427 

|zl)-1.25 

431  SUM 

-429 

408 

iz2) - .01 

432  ORR 

430 

431 

outside  outer  rectangle? 

433  INV 

432 

true  if  Inside  outer  rectangle 

434  INV 

422 

true  if  outside  zero  rectangle 

1.0000 

435  SRS 

4  34 

433 

Identify  and  loop  tliru  regions 

436  AND 

4  35 

434 

compare  loop  and  location 

437  IFE 

436 

421 

409 

Uz~0  until  outside 

0.0000 

439  CON 

INPUT  (c/J) 

440  HUL 

4  39 

450 

(c/J)Tz 

0  01818/.3 

441  CON 

INPUT  (1/)) 

0 . 0000 

443  CON 

INPUT  (f/J) 

444  riUL 

-443 

450 

(-f/j)Tz 

445  SUM 

72 

244 

437 

sum  of  control  torques 

446  HUL 

14 

413 

404 

4LUcz 

447  MUL 

18 

404 

2Ucz 

448  SUM 

6 

446 

1 4  4  LUc  z 

449  DIV 

448 

44  7 

(l44LUcz)/(2Ucz) 

U  0000 

450  CON 

INPUT  Tz 

Fi'j’ure  30.  TUTSI M  Simulation  Program  (Third  Page) 
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APPENDIX  B.  SIMULATION  PLOTS 


Figure  31.  Siimilalioii  Plots  for  Baseline  No  Target 


Figure  32.  Baseline  No  Target 
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Figure  36.  Baseline  Case  IB 
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Fi"uie  38.  Baseline  Case  1C 
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Figure  39.  Baseline  Case  ID 
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Fi<’iiie  40.  Baseline  Case  ID 
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Figure  44.  Proposed  No  Target 
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Figure  56.  Proposed  Case  ID 
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Figure  61.  Proposed  Case  IB  End  of  Mission 
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Figure  62.  Proposed  Case  IB  End  of  Mission 
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B 1 oc  kNo , 

Plot-HINl mum . 

Plot  -NAXlmuni; 

Comniriit 

201  . 

-0.0870000  . 

0.0870000  : 

INPUT  yl 

INITIAL 

CONDITION 

200  . 

-0.0870000  , 

0.0870000  : 

INPUT  y2 

INITIAL 

CONDITION 

yi:lHrUT  yZ  IHITIfIL  COMP  IT  ION 
0.0070000  '  ^ 

O.0bOGO0O 

0.0522000 

0.031UOO0 

0.0171000 

0.0000 

-0.0171000 

-0.0310000 

-0.0522000 

-0.063GOOO 

-0.0070000  - 1 _ L. 


CftSElC.l  EHD  PITOI  mOPOSED 


_ I _ L_ _ < _ I - U - 5 

-0.0O7OO0INrUT  ul  INITIAL  CONDITION 


0.0070000 


Block  Ho,  Plot  KIN  Innim 
GOl  ,  -0.0870000 

GOO  ,  -0,1750000 


r I o t ■ MAX  I  mum 
0.0870000 
0.  1750000 


ComiBPnt 

INPUT  zl  INITIAL  CONDITION 
INPUT  z2  INITIAL  CONDITION 


Yl:  INPUT  7.Z  INITIAL  CONDITION 
0 . 1750000  [  '  ^ 


0.1100000 
0 . 1050000 
0.0700000 
0.0350000 
0.0000 
-0.0350000 
-0.0700000 
-0.1050000 
-0.1100008 
-0.1750000 


CASElC.l  END  VAU  rnOPOSED 


-0.007000 INPUT  zl  INITIAL  CONDITION 


0.0070000 


Fi”iii  e  6-}.  Proposed  Ca.sc  1C.  I  End  of  Mission 


Figure  65.  Piuposcd  (  ase  ID  Ciui  of  Mission 


Fimirc  6(').  Proposed  (  ase  ID  End  of  Mission 


ri”iire  67.  Proposed  Case  2  End  of  Mission 


Yl:IHrUT  yZ  INITIAL  CONDITION 
0.0070000  ' 

0.O&9GOO0 

0.0S2Z0O0 

0.03-10000 

0.017-1000 

0.0000 

-0.017-1000 

-0.03-10000  - 

-0.0S22000 

-O.0A960OO 

-0.0070000  - 1 _ L 


C0SE2  END  rnai  rnoroscD 


-0.007000 INPUT  yl  INITIAL  CONDITION 


0.0070000 


B 1 ockNo , 

rio  t  •  M I  Nl  tmim , 

Plot  -MAXIimini; 

Comment 

/<01  , 

-0.0870000  . 

0.0870000  : 

INPUT 

zl 

INITIAL 

COimiTION 

GOO  , 

-0.0870000  , 

0.0870000  ; 

INPUT 

z2 

INITIAL 

CONDITION 

VI : INPUT  t2  INITIAL  CONDITION 
0.0070000  ^ 

0.0O9G000 

0.0522000 

0.0310000 

0.0171000 

0.0000 

-0.0171000  - 

-0.0310000 
-0.0522000 
-0.069G000 

-0.0070000  - 1 _ L 


CASE2  END  YOU  rOOPOSED 


-0.0O700OINPUT  zl  INITIAL  CONDITION 


0.0070000 


Fi”ui  e  68.  Proposed  Case  2  End  of  Mission 


LIST  OF  REFERENCES 


1.  Best  and  Final  Offer;  Responses  to  Questions  and  RFP  Amendment  7  Book  2, 
McDonnell  Douglas  Astronautics  Company-Space  Station  Division,  25  September 
19S7. 


2.  Space  Siaiion  Projects  Office  ll'ork  Package  2  Request  for  Proposal 
^9-BF2-51-7-l P-POl j ,  National  Aeronautics  and  Space  Administration,  Johnson 
Space  Center.  May  19S7. 

3.  Hamilton.  Sir  W.R..  Elements  of  Quaternions,  pp.113-114.  London:  Longmans. 
Green  and  Co..  1S66. 

4.  Wertz.  James  R.,  Spacecraft  Attitude  Determination  and  Control,  Holland:  D. 
Reidel  Publishing  Company.  1986. 

5.  Halliday.  D.  and  Resnick.  R..  Fundamentals  of  Physics  3rd  ed.,  pp.  236-24(').  Canada: 
John  \\'iley  cC  Sons.  Inc..  1988. 

6.  .Martin  Marietta  .Aerospace  Report  M.MU-SE-17-46.REV  A.,  Manned  Manuevering 
Unit,  Space  Shuttle  Program  User's  Guide,  pp.  54-59.  July  1985. 

7.  Goldstein.  H..  Classical  Mechanics.  2nd  ed..  pp.  197-198,  Reading.  Vlassachusetts: 
Addison-^^'esIey,  1981. 

8.  Kirk.  D.  E.,  Optimal  Control  Theory:  An  Introduction,  Prentice-Hall  Inc., 
Englewood  ClilEs,  New  Jersey,  1970. 

9.  Reynolds.  W.E..  and  Wolf.  J.,  TUTSIM  Users  Manual,  Tutsim  Products.  Palo  .Alto. 
California.  September  1988. 


100 


INITIAL  DISTRIBUTION  LIST 

No.  Copies 

1.  Defense  Technical  Information  Center  2 

Cameron  Station 
Alexandria,  VA  22304-6145 

2.  Library,  Code  0142 
Naval  Postaraduate  School 
Monterey,  CA  93943-5002 

3.  Commander 
United  States  Space  Command 
Attn:  Technical  Library 
Peterson  ALB.  Colorado  80914 

4.  Commander 
Naval  Space  Conunand 
Attn:  Code  N3 
Dahlgren,  Virginia  22448 

5.  Director.  Navy  Space  Systems  Division 
Chief  of  Naval  Operations 
(OP-943) 

Washington,  DC  20350-2000 

6.  Superintendent  1 

Attn:  Chairman.  Code  62 

Department  of  Electrical  and  Computer  Engineering 
Naval  Postgraduate  School 
Monterey,  California  93943 

7.  Superintendent  4 

Attn:  Professor  R.  Panholzer,  Code  72 

Space  Systems  Academic  Group 
Naval  Postgraduate  School 
.Monterey,  California  93943 

8.  Superintendent  2 

.Attn:  Prof  Jelf  B.  Burl.  Code  62B1 

Department  of  Electrical  and  Computer  Engineering 
Naval  Postgraduate  School 
.Monterey,  California  93943 

9.  Superintendent  1 

Attn:  Prof  11.  A.  Titus.  Code  62Ts 

Department  of  Electrical  and  Computer  Engineering 
Naval  Postgraduate  School 
.Monterey,  California  93943 


2 


2 


") 


101 


1 


10.  Superintendent 
Attn:  Dr.  George  J.  Thaler.  Code  62Tr 
Department  of  Electrical  and  Computer  Engineering 
Naval  Postgraduate  School 
Monterey.  California  93943 

11.  McDonnell  Douglas  Space  Systems  Company  1 

Attn:  Mr.  John  Harduvel 

A95-J845;  M  S  11-3 
5301  Bolsa  .Avenue 

Huntington  Beach,  California  92647-2048 

12.  .Mr  C.  E.  Whitsett  1 

N.AS.A.  Johnson  Space  Center 

EC5 

Houston,  Texas  77058 

13.  Commander  1 

Naval  Air  Force.  U.S.  Atlantic  Fleet 

Attn:  CAPT  T.C.  Eackey.  Code  03 
Norfolk.  Virginia  2351  1 

14.  Commander  1 

VM  FAT- 101 

.Attn:  ET  JelTrcy  W.  Hansen 
MC.AS  El  Toro.  California  92630 

15.  Commander  2 

Helicopter  Combat  Support  Squadron  Six,  LP4 

.-\ttn:  l.CDR  Daniel  L.  llansen 
N.AS  Norfolk.  Virginia  2351  1 


102 


